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We provide a general quantum theory to describe the couphng of hght with the motion of a 
dielectric object inside a high finesse optical cavity. In particular, we derive the total Hamiltonian 
of the system as well as a master equation describing the state of the center of mass mode of the 
dielectric and the cavity field mode. In addition, a quantum theory of elasticity is used in order 
to study the coupling of the center of mass motion with internal vibrational excitations of the 
dielectric. This general theory is applied to the recent proposal of using an optically levitating 
nanodielectric as a cavity optomechanical system P, 0]. On this basis, we also design a light- 
mechanics interface to prepare non- Gaussian states of the mechanical motion, such as quantum 
superpositions of Fock states. Finally, we introduce a direct mechanical tomography scheme to 
probe these genuine quantum states by time of flight experiments. 



I. INTRODUCTION 

The field of optical trapping and manipulation of small 
neutral particles using the radiation pressure force of 
lasers was originated in 1970 by the seminal experiments 
of Ashkin [3]. Over the course of the next 40 years, the 
techniques of optical trapping and manipulation have 
stimulated revolutionary developments in the fields of 
atomic physics, biological sciences, and chemistry [4]. In 
physics, the progress in optical cooling and manipulation 
of single atoms opened up a plethora of novel perspec- 
tives. The precise control over the atomic degrees of free- 
dom has created applications ranging from atom inter- 
ferometry [5], quantum simulations of condensed matter 
systems with ultracold gases [6] , and the implementation 
of quantum gates for quantum computation purposes [t^ . 

More recently, the possibility to apply the techniques of 
optical cooling and manipulation to the mechanical de- 
gree of freedom of larger objects, such as micromirrors 
or cantilevers, has established a very active research field 
- cavity quantum optomechanics jsl-ll^ . Future applica- 
tions range from ultra-high sensitivity detectors of mass- 
or force [15|, [l6| and quantum transducers for quantum 
computation purposes p!7l-[2Q|, to their potential of be- 
ing an ideal testbed for the investigation of fundamental 
aspects of quantum mechanics, such as the quantum- to- 
classical transition [2l|, [22| . In most optomechanical sys- 
tems, the mechanical oscillator is unavoidably attached 
to its suspension providing a thermal contact that limits 
the isolation of the mechanical motion - thus prevent- 
ing longer coherence times. A potential improvement to 
better isolate the system is to use optically levitating 
nanodielectrics as a cavity quantum optomechanical sys- 
tem Q (see also [23|, uM)- This consists in optically 
trapping a nanodielectric by means of optical tweezers in- 
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FIG. 1: Schematic representation of the setup. A nanodielec- 
tric is confined by optical tweezers which provides a trapping 
frequency of cut- The nanodielectric is placed inside an opti- 
cal cavity with resonance frequency cjc, decay rate /€, and is 
driven by a laser at a frequency uol- 



side a high finesse optical cavity, see Fig. [T] for an illustra- 
tion. Using standard optomechanical techniques p5l-[29j. 
the center of mass (CM) motion could be cooled to its 
quantum ground state in the harmonic potential created 
by the optical tweezers. Owing to the fact that it is lev- 
itating, the dielectric is not attached to any mechanical 
object, allowing for a very good thermal isolation, even 
at room temperature. More recently, both theoretical 
and experimental research along this direction has been 
reported. In [l6|, levitating nanospheres placed close to 
a surface have been proposed to test forces at very small 
scales in order to explore corrections to the Newtonian 
force. Remarkably, an experiment measuring the instan- 
taneous velocity of the Brownian motion of a particle, 
a glass bead levitated in air, has been reported in (soj . 
Besides, other aspects have been investigated, such as 
the possibility of Doppler cooling a microsphere [3l|, a 
scheme to measure the impact of air molecules into the 
nanodielectric [32|, and the possibility to use a ring cavity 
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to cool and trap polarizable particles [33| . It can thus be 
foreseen that a new generation of exciting experiments, 
aiming at bringing levitating dielectrics into the quan- 
tum regime, will eventually take place in the near future. 
Indeed, from a broad perspective, this project aims at ex- 
tending the techniques developed during the last decades 
of optical cooling and manipulation of atoms (e^. like 
in cavity QED with single atoms and molecules [3j437|), 
back to the nanodielectrics that were first used in the 
times of birth of optical trapping (ssMioj . This experi- 
mental challenge, if successful, would allow to test quan- 
tum mechanics at unprecedented scales. On this basis, a 
general quantum theory to describe and predict the phe- 
nomena to be encountered in these potential experiments 
is timely. 

In this article, we aim at contributing to this goal by de- 
veloping a general quantum theory to describe the cou- 
pling of light with the mechanical motion of dielectrics 
in high finesse optical cavities. Starting from the total 
Hamiltonian of the system, we derive a master equation 
which takes into account the effects caused by the scat- 
tering of light: decoherence of the mechanical motion, 
decrease of the cavity finesse, and a non- negligible renor- 
malization of the scattering force. Additionally, we uti- 
lize a quantum elasticity theory to describe the effect of 
elastic deformations of the dielectric object. This pro- 
vides us with a quantitative expression of the coupling 
between the center of mass motion and the internal vi- 
brational modes. This theory is applied to the particular 
proposal of cavity optomechanics with levitating nanodi- 
electrics [1, 2]. We then focus on "post ground state" [80| 
optomechanics, and design a light-mechanics interface to 
prepare non-Gaussian states of the mechanical system, 
such as a quantum superpostion of Fock states. In par- 
ticular, we develop three different protocols with different 
features, together with a formalism which is required to 
describe these input-ouptut problems in the Schrodinger 
picture. The non- Gaussian light-mechanics interface can 
be interpreted as an effective way to have non-linear ef- 
fects in optomechanical systems [4l|. Finally, we in- 
troduce a scheme to perform direct full tomography of 
the mechanical state by imaging the nanodielectric after 
some time of flight. 

The article is organized as follows: in Section [iTl we pro- 
vide a detailed summary of the results. The derivation of 
the total Hamiltonian as well as the master equation of 
the theory is addressed in Sec. IIIII (some details regarding 
light scattering are left to App. [Aj). The derivation of the 
trapping using optical tweezers and the optomechanical 
coupling can be found in App. [B1 The theory part is 
completed by the introduction of the quantum theory of 
elasticity in Sec.[lVl We briefly discuss ground state cool- 
ing in App. \C\ and the typical experimental parameters 
in App. [El The second part of the article, the description 
of the protocols, can be found in Sec. El (and in App. [P]). 
where three different ways to interface light with the cen- 
ter of mass motion of the nanodielectric are introduced. 



The article is rounded off by a proposal to perform full 
tomography of the mechanical state in Sec. Wl\ and we 
finish by stating the conclusions and discussing further 
directions in Sec. IVIIi 



II. SUMMARY OF RESULTS 



This section aims at providing a general roadmap, sum- 
marizing the results presented in this manuscript without 
providing the proofs or mathematical derivations. 



A. Theory 

Here we develop a quantum theory to describe the cou- 
pling of light with the mechanical motion of dielectric 
objects in high finesse optical cavities, see Fig. [H The 
assumptions that are made in the theory are the follow- 
ing: 

1. The dielectric object has a constant relative dielec- 
tric constant e^, as well as, a homogeneous density 
P- 

2. The electric fields are assumed to be scalar, that 
is, we do not consider polarizations. We assume a 
three dimensional field for the free electric field out- 
side the cavity, and a one-dimensional field along 
the cavity axis for the output field of the cavity. 

3. The object is assumed to be absorption- free and 
therefore only elastic scattering processes are con- 
sidered. The effects of light absorption are thus 
neglected [l|,|2[. 

These assumptions are made in order to ease the deriva- 
tion of the theory and do not imply any fundamental 
simplification. Indeed, non-homogeneity and polariza- 
tions can be incorporated easily. Moreover, as it is shown 
in the quantum theory of elasticity, the center of mass 
mode is decoupled from the internal vibrations for suffi- 
ciently small objects, and therefore can be treated inde- 
pendently. 

Thus, the total Hamiltonian of the system can be written 
as a sum of free (interacting) terms, labeled with the 
super-index f (i), 

^tot = + ^cav + ^out + ^free + ^cav-out + ^diel' (1) 

The first term = /2M is the kinetic energy of the 
center of mass position along the cavity axis. The mo- 
tion along the transverse direction of the cavity is not 
relevant in the theory jsij. The energy of the cavity 
mode a is given by i^^av — ^cO^^cl (we assume h = 1), 
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where uOc is its resonance frequency. The energy of 
the free modes is given by i^l^ee ~ / <^k|k|a'''(k)a(k), 
and the energy of the ouptut modes of the cavity by 
^out — dujuja\{(jj)aQ{uj) [82|. The interaction between 
the cavity mode and the output modes is described by 
the usual term i^cav-out = i /o^ duj'^iuj) [o^aQiu) — H.c), 
where the couphng strength is approximated by ^{uo) ^ 
a^/tt, around the resonance frequency, where is the de- 
cay rate of the cavity (42|. 

The last term i^^^gi is the crucial one describing the inter- 
action between the electric field and the dielectric object. 
In the most general form, it can be written as 



^diel 



y(r) 



dxP(x)£:(x), 



(2) 



where P(x) is the polarization of the object and the in- 
tegration is performed over the volume of the dielectric 
object V with center of mass coordinate r. This Hamil- 
tonian is the starting point for the theoretical discussion. 
Assuming P(x) = ap£^(x), one obtains 



^diel — 
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y(r) 



(3) 



where Cc = 3(e^ — l)/(e^ + 2), being the relative dielec- 
tric constant of the nanodielectric. This can be obtained 
by connecting the quantum expression of the polarization 
field in the object with the classical relation. This part of 
the Hamiltonian is the key ingredient of the theory, and 
applies for any size and shape of the object. The total 
electric field inside the object can be now written as a 
sum of three parts E((x) = £^cav(x) + £^free(x) + ftw(x), 
where ^cav(x) contains the cavity modes, £^free(x) the 
free modes, and ftw(x) is the classical part of the electric 
field describing the optical tweezers. By plugging £^(x) 
into Eq. ([3]) the following terms are obtained: 

1. [^tw(x)]^ creates a harmonic trap with a frequency 
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where / is the laser intensity, p the density of the di- 
electric object, c the speed of light, and Wt the laser 
waist. For the typical experimental parameters dis- 
cussed in App. [E] it is of the order of MHz. This 
field provides the harmonic trap of the mechanical 
oscillator with the Hamiltonian outb^b, where b (b^) 
is the annihilation (creation) operator of the center 
of mass phonon mode along the cavity axis. 

2. The cavity field [£^cav(x)]^ gives rise to the optome- 
chanical coupling goa^a{b^ +6), where the coupling 
strength is given by 
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Here zq = {2Mujt)~^^'^ is the zero point motion of 
the ground state, Vc = LiiW'^ /A the cavity volume, 
L the cavity length, and Wc the laser waist of the 
cavity. For the experimental parameters discussed 
in App. [El qq is of the order of tens of Hz. 

3. The contribution 2£^free(x)(£^cav(x) +ftw(x)) is re- 
sponsible for scattering processes. This term de- 
scribes the process of elastic scattering of cav- 
ity photons, as well as photons of the tweezers, 
into free modes. The term 2^tw(x)£^cav(x) leads 
to scattering events already taken into account in 
2£^free(x)ftw(x) as wcll as a shift in both the trap- 
ping frequency and the equilibrium position of the 
object. 

4. The term [E'free(x)]^ yields a negligible coupling be- 
tween the center of mass mode and the vacuum 
fluctuations of the electromagnetic field which is 
negligible. 

Starting from the total Hamiltonian i^tot and the terms 
given by the total electric field £^(x), one can derive a 
master equation describing the state of the cavity mode a 
and the mechanical mode 6, given by the density matrix 
p, by tracing out the free modes a(k) and the ouptut 
modes ao(cj). The master equation is given by: 

with the following contributions: 

1. The optomechanical Hamiltonian in the non- 
displaced frame [83| is given by 



i^OM = ojtb'^b + uJcCi^ a + gQa^a{b'^ + 6), 



(7) 



(5) 



describing the coherent coupling between the cavity 
mode and the mechanical mode. 

2. The dissipation term 

>^cav[p] = i^{2apa^ — ap — pa^a)^ (8) 

describing the photon losses, at a rate z^, due to the 
imperfection of the cavity mirrors. 

3. The new cavity field dissipation term 

^sc[p\ = t^sc(2apa^ — a^ap — pa^a), (9) 

due to losses, at a rate tv^c caused by the scattering 
of cavity photons out of the cavity. Although the 
theory is valid for any size of the object, we pro- 
vide in the article the expression of hCsc for objects 
smaller than the optical wavelength. Indeed, in or- 
der to keep the high finesse of the cavity, that is, 
^sc/^ < the objects have to be of the order of 
100 nm in case of spherical objects [2[. 
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4. The mechanical diffusion term 

Vm[p]=rsc[b + b^,[b + b^,p]], (10) 

which, although it does not yield mechanical damp- 
ing, does generate decoherence of the motional me- 
chanical state due to light scattering. We also pro- 
vide in this article the expression of Fgc for sub- 
wavelength objects, which contains the contribu- 
tion of both the optical tweezers and the cavity 
field. For spherical objects of the order of 100 nm, 
l/Fsc ~ 0.1 ms. 

5. Finally, an additional coherent term Hm is ob- 
tained, which renormalizes the optomechanical 
Hamiltonian due to QED effects. The effects of 
this term are discussed in more detail in j43} . 

This theory is complemented by a quantum theory of 
elasticity. Starting from the classical expression of the 
Fagrangian density of an elastic object, the deformable 
field is expressed in terms of normal modes for the case of 
a vanishing external potential. Then, by plugging in the 
external potential given by the light matter interaction, 
one obtains an expression describing coupling between 
the normal modes. This can be quantized canonically 
and provides a quantum description of the coupling of 
the center of mass mode with the internal vibrational 
modes, as well as the coupling of the light with the inter- 
nal modes. This theory can be applied to objects at the 
micron scale, where the internal modes have frequencies 
of the order of 10^^ Hz, and thus are decoupled from the 
10^ Hz center of mass mode. This allows us to adiabat- 
ically eliminate the internal modes, merely leading to a 
renormalization of the trapping frequency. This correc- 
tion is many orders of magnitude smaller than ujt and 
consequently represents a negligible effect. This justifies 
the separate treatment of the center of mass degree of 
freedom in the developed theory which is applicable to 
objects at the the micron scale. 



B. Protocols 

In the second part of the article we focus on how to bring 
these systems into the quantum regime. In particular we 
design a light-mechanics interface which consists in in- 
jecting non-Gaussian states of light, such as superposi- 
tion of Fock states, to the mechanical oscillator. 

First of all, as usually done in optomechanics, we de- 
scribe the effect of driving the cavity with a strong driv- 
ing field at frequency ujl [26, 28, 44]. We transform the 
total Hamiltonian of the system into a displaced frame 
which describes the states on top of the steady state of 
the cavity field, the mechanical state, as well as the out- 
put modes of the cavity. We particularly emphasize the 
displacement that one needs to do in the output modes 



in order to be able to describe input-output problems in 
the Schrodinger picture. The main change in the Hamil- 
tonian is in the optomechanical coupling. In particular, 
when the driving is red-detuned A = Uc — ujl = ^t, the 
coupling term, in the resolved sideband regime and after 
the rotating wave approximation, has the beam-splitter 
interaction form g{a^b-\- ab^). Here, g = ^/n^go is an ef- 
fective optomechanical coupling, which is enhanced by 
the square root of steady state cavity photons. This 
allows one to reach the strong coupling regime g > n. 
When driving the cavity with the blue-detuned field 
/S. = uJc — ojl = —^ti one induces the two mode squeezing 
interaction term g{a^})^ ^ah). With these tools, we design 
and derive different protocols to perpare non- Gaussian 
states. 

The first protocol, called the refiected one-photon, con- 
sists in sending one resonant photon on top of the driving 
field and measuring the refiected part. More specifically, 
the cavity is driven with a red-detuned field in order to 
induce the beam-splitter interaction. The mechanical ob- 
ject is assumed to be in its ground state. Now, on top 
of the driving field, a one-photon pulse centered at the 
resonance frequency is sent into the cavity. Impinging 
the cavity, part of it enters and part is refiected. At the 
time t/i, where the part of the beam that has entered the 
cavity is transferred to the mechanical oscillator through 
the beam-splitter interaction, the light field is switched 
off. Consequently, the light mode corresponding to the 
refiected photon is entangled with the mechanical sys- 
tem inside the cavity. We can obtain the exact form 
of the state by solving the input-output problem in the 
Schrodinger picture. The state in the displaced frame is 
given by 

/oo 
c{u:,th)al{u:)<h:\QQQ), 

(11) 

where In^na^^) describes a state with phonons, Ua pho- 
tons, and all the output modes in the vacuum state. Here, 
the coefficients ci){t) and c(cj, t) are obtained analytically. 
This makes clear that, by measuring the quadrature of 
the output mode of the photon, one prepares a super- 
position state of zero and one phonon with coefficients 
given by the outcome of the measurement. Some techni- 
cal issues are addressed in detail for this protocol in the 
manuscript, such as the fact that in the original frame, 
the state \il){th)) is displaced by a considerable amount. 
This makes it challenging to obtain a significant signal- 
to-noise ratio in the measurement of the output mode. 

An extension of the refiected one-photon protocol is the 
perfect mapping protocol. In this protocol, the possibil- 
ity to time-modulate the laser intensity, and consequently 
the optomechanical coupling, is exploited. Then, by im- 
posing the condition that the output field, with the trans- 
formed Hamiltonian, is zero, we can obtain the equation 
of motion for the optomechanical coupling g(t). Its solu- 
tion yields a modulation of g{t) such that the light pulse 
sent on top of the driving field is perfectly absorbed and 
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therefore the non- Gaussian state of the hght is trans- 
ferred to the mechanical oscihator. In this section we 
also discuss some technical details regarding the trans- 
formation of the Hamiltonian that has to be performed 
carefully since there are time-dependent displacements. 

The two protocols require a moderately strong coupling 
g ^ tv. As an alternative, we also derive a protocol, called 
teleportation in the bad-cavity limit, which does not re- 
quire the strong coupling regime. Once the mechanical 
oscillator is in the ground state, it consists in driving the 
cavity with a blue-detuned field, such that the two mode 
squeezing interaction is induced inside the cavity. This 
Hamiltonian creates a two mode squeezed state between 
the mechanical oscillator and the light field leaking out 
of the cavity. The squeezing parameter of this state is a 
measure of the degree of entanglement. This entangled 
state can then be used to teleport a non-Gaussian light 
state from outside of the cavity to the mechanical oscilla- 
tor. In this section, we also discuss in detail the effect of 
the driving field and the possibility to choose the initial 
state to be teleported in order to prepare a particular 
state in the mechanical oscillator. 

This part of the article is concluded by providing a di- 
rect method to perform full tomography of the state of 
the mechanical oscillator. In general optomechanical sys- 
tems, tomography can in principle be done by coupling 
the mechanical resonator to a well-controlled quantum 
system {e.g. a qubit), and then measuring the quantum 
system. The method we provide here performs direct to- 
mography of the mechanical oscillator. It is well known 
that by measuring the rotated quadrature phase operator 

X{0) = e'%^ ^e-'^b, (12) 

for all 6>, one can reconstruct the Wigner function and 
therefore obtain all the information about the state of the 
harmonic oscillator [45^. Our protocol achieves that by 
measuring the position of the nano-dielectric after some 
time of fiight. In the Heisenberg picture, the momentum 
operator in the harmonic potential evolves like 

p{t) = ip^ih^e'^^' - &e--'*), (13) 

where Prn = {Mujt/'^Y^'^ • Thus the momentum operator 
p{t)/Pm = ^{^tte + '7r/2) is directly related to the ro- 
tated quadrature phase operator. By switching off the 
optical tweezers at tg, letting the nanodielectric fall, and 
measuring the position at some later time t/, one obtains 
z{te-\-tf) ~ {tf —te)p{te)/M^ which is a measurement of 
the momentum operator. By repeating the experiment 
at different times te full tomography of the mechanical 
state can be performed. In this section, we discuss the 
experimental parameters and an extension of the proto- 
col to amplify the oscillation prior to the time of fiight. 



III. TOTAL HAMILTONIAN AND MASTER 
EQUATION OF THE THEORY 

In this section we develop the main part of the theory to 
describe the coupling of light with the mechanical mo- 
tion of dielectric objects in high finesse optical cavities. 
In particular, in Sec. IIII Al we derive the total Hamilto- 
nian. We focus on the light-matter interaction term in 
Sec. IIII B\ and make the connection between the micro- 
scopic theory with the macroscopic parameters such as 
the dielectric constant of the object. In Sec. IIII CI we 
show how to obtain the optomechanical Hamiltonian for 
the case of levitating objects inside a cavity. Finally in 
Sec. HHP I we derive the effects of scattering of light em- 
bedded in the total Hamiltonian of the theory. Indeed, 
we show how to derive a master equation in order to 
describe the time evolution of the state of both the cav- 
ity and mechanical mode. We provide the quantitative 
expression of the light scattering parameters for subwave- 
length spheres. 



A. Setup and general Hamiltonian 

We consider a dielectric object with center of mass posi- 
tion r and a dielectric constant 

where V(r) is the spatial region of the object, centered 
at r, with volume density p and mass M = pV. The 
homogeneity of the dielectric constant inside the object 
is chosen for simplicity, the non-homogeneous case can 
be incorporated easily. As shown in Sec. HVl the center 
of mass degree of freedom of dielectrics at the micron- 
scale is decoupled from its relative modes. Hence, we 
will only consider the center of mass degree of freedom 
in the following analysis. We suppose that the dielectric 
object is inside an optical cavity. We define the cavity 
mode, characterized by an annihilation (creation) opera- 
tor a {a^), the modes coupled to the cavity mirror ao{uj) 
{al{uj)), we call them output modes, and the other free 
modes of the electromagnetic field a(k) (a^'k)). We use 
a ID theory to describe the modes coupled to the cavity 
mode, therefore denoting them hy uj = k (we use c = 1 
in the protocols part of the article) and a scalar 3D the- 
ory for the rest of the modes. The effects of polarization 
can also be easily incorporated but will be neglected for 
simplicity. 

The total Hamiltonian of the dielectric object inside the 
optical cavity can be written as 

^tot =^m + ^cav + ^out + ^free + ^cav-out + ^diel' 

(15) 

The superscript / (i) labels free (interacting) terms. The 
first term, = /2M describes the kinetic energy 
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of the center of mass mode along the cavity axis. The 
motion along other directions is not considered since it 
decouples from the motion along the cavity axis, as it 
will become clear in the following. Note however that 
in analogy to trapping and cooling of ions, these other 
modes are assumed to be cooled by external means {e.g. 
by feedback coohng) in order to make the trap stable (see 
a recent article for a 3D ground state cooling scheme 
based on using different cavity modes [32|). The next 
three terms describe the free radiation parts 

^cav = ^ctt'^a, (16) 

^out = / duju;al{uj)ao{oj), (17) 
Jo 

Hi,, = I ^k|k|at(k)a(k), (18) 

of the cavity mode a with the mode frequency cjc, the out- 
put modes ao{uj) and the free modes a (k). As it is usually 
done in the context of cavity QED [46| , we double-count 
some of the modes by considering the ones coupled to 
the cavity mode separately. However, since they have 
zero measure, this does not affect the correctness of the 
description. The interaction between the cavity mode 
and the free modes is described by (42^ 

^cav-out = i / cl^lH (a^aoiuj) - H.c.) , (19) 
Jo 

where the coupling strength is described by j{oj) and can 
be assumed to be constant over a large frequency interval 
centered around lJc with a value ^(co) = ^J~KpK^ where k 
is the decay rate of the cavity [42]. 

Finally, describes the interaction between the light 

field and the dielectric object which can be written as 

^iiei = d^P{^)E{^). (20) 

Here, P(x) is the polarization field and the volume inte- 
gral is performed over the volume of the object V around 
the center of mass position r. This term is the central 
equation in the rest of subsections: in Sec. IHI B\ we de- 
velop this interaction term by relating the polarization 
field with the electric field; in Sec. IHI C\ we consider the 
proposal of using an optically levitating nanodielectric 
as a quantum optomomechanical system, and derive the 
optomechanical Hamiltonian; and in Sec. IHID^ we show 
how this term can be used to derive the effects induced 
by light scattering. 

As a final remark and for later convenience, let us define 
the light mechanics (LM) and light-cavity (LC) part of 
the Hamiltonian as 

^LM = ^cav + ^free + ^dieb ^^^^ 
HlC = ^out + ^cav-out' 

such that i^tot = i^^ + ^lm + i^LC- 



B. Light-matter interaction Hamiltonian 

Let us here focus on the key part of the total Hamilto- 
nian, the light-matter interaction term i^^i^i, Eq. ([2Q|) . 
First of all, note that for the typical light intensities, the 
polarization field responses linearly to the electric field, 
such that P(x) = apE{x.). The parameter ap can in 
principle be computed by performing a quantum theory 
of the object by considering its atomic properties. How- 
ever, this involved task is not necessary since one can 
relate ap to macroscopic properties of the object, such as 
the dielectric constant e^. Comparing the resulting re- 
lation between the polarization and the electric field for 
the macroscopic [43, |48| and microscopic case: 

macroscopic: P(x) = Scq— -£^(x) = eceo£^(x) 

Cf -\- 2i (22) 

microscopic: P(x) = apE{x.) 

one can identify the microscopic constant to the macro- 
scopic one, 

c^p = ecCo. (23) 

Then, plugging this back into the Hamiltonian Eq. ([2Q|) . 
we obtain the final form of the light-matter interaction 
Hamiltonian 

HL^ = / rfx[E(x)]2. (24) 

In App. [Al we show how from this Hamiltonian one can 
derive the scattering equation that can be used to com- 
pute the electric field inside the object. 



C. Optomechanical Hamiltonian 

The expression for i^^^gi obtained in the above section, 
see Eq. ([24|) . is the key ingredient to develop our the- 
ory. Let us now apply it to the particular proposal of 
using levitated objects in a cavity as an optomechanical 
system [l|, Q . This experimental setup consists of an ex- 
ternal optical tweezers as well as a laser driving the cavity 
at frequency cjl, see Fig.[TJ The total electric field inside 
the object can be decomposed into 

^(x) = ^cav(x) + ^free(x), (25) 

where the £^cav(x) is the cavity electric field and £^free(x) 
the free electric field. The free electric field contains a 
classical part due to the optical tweezers generated by the 
laser, which can be incorporated as £^free(x) Efree (x) + 
ftw(x), where ftw(x) describes the optical tweezers |49j , 
see App. [B] for its expression. The implementation of the 
driving laser is carefully done in Sec. IV Al where, in order 
to keep the structure of the total Hamiltonian, we will 
have to displace all the output modes ao{uj) as well as 
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the cavity mode and the mechanical mode. Let us here 
discuss the terms that wih be obtained when plugging 
the total electric field 



£;(x) = £;cav(x) + £^free(x) + ^tw(x) 



(26) 



in i^^iei' (EH)- -^y doing so, one obtains six differ- 
ent terms, which have been discussed in the summary of 
results section, and therefore are only summarized here. 
The term [ftw(x)]^ which accounts for the harmonic trap- 
ping, see App. [HI with a trapping frequency 



4ec_/_ 

pc wr 



(27) 



where / is the field intensity, Wt ^ XjiiiM) the laser 
waist, M the numerical aperture, and k the wave vector 
number. This allows us to quantize the CM motion along 
the z axis as z = z^(\)^ + h\ where = (2MijJi)~^l^ . 
The term [E'cav(x)]^ describes the coupling of the cavity 
mode and the motional state, see App.lBl By considering 
the center of mass position of the object to be placed 
at the maximum slope of the standing wave one obtains 
the standard optomechanical coupling g^a^a{\)^ +6). The 
coupling strength is given by 



^0 



V 



(28) 



where V is the volume of the object, Vc = LttW^/A 
the cavity volume, L the cavity length, and Wc = 
[AI//(27r)]^/^ the waist of the cavity field. This term also 
yields a shift of the resonance frequency of the cavity, 
see App.lBl The two scattering terms 2£^free(x)(ftw(x) + 
£^cav(x)), which describe the scattering of cavity pho- 
tons and the laser light from the optical tweezers, are 
addressed in Sec. IIIIDl The term 2£^cav(x)ftw(x), which 
yields a shift of both the trapping frequency and the equi- 
librium position as well as some scattering processes al- 
ready taken into account in the term 2£^free(x)ftw(x), is 
discussed in the next section. Finally, the term [£'free(x)]^ 
accounts for a negligible coupling of the center of mass 
motion with vacuum fiuctuations of the electromagnetic 
field. 



Hence, the total Hamiltonian can be now written as 



Hqm + Hsc + i^LC + 



sh 7 



(29) 



where Hqm is the standard optomechanical Hamiltonian 
^OM = ojtb^b + ujcd^a + goa^a{b^ + b). (30) 



The term i^Lc was already introduced in Eq. ()2T|) and 
we have defined the scattering Hamiltonian 

Hsc = i^free - ^cCO / C^x£;free (x) (^^cav (x) + ftw(x)), 
Jvir) 

(31) 

which is studied in Sec. IHIDi The shift term is given 
by i^sh = -ecCo / £^cav(x)ftw(x)c^x. Finally, note that 



by tracing out the output modes of the cavity, ao{uj) in 
i^LC? one obtains the usual master equation 



p{t) = i[p, HoM + Hi^c + i^sh] + Cc^Ap] (32) 



where 



^cav[p] = i^{2apa^ — a^ap — po) a) 



(33) 



This term describes the photon losses (at the decay rate 
of the cavity /^) through the end mirrors of the cavity, 
and is treated in detail in, e.g., (42j . 



D. Light Scattering 



In the present setup both trapping and cooling is 
achieved by light, yielding an optomechanical system 
without thermal contact to other mechanical objects. 
However, the effect of scattering of light has to be consid- 
ered [2] . In this section we provide a framework to study 
elastic light scattering within a quantum theory. While 
the derivation of this framework for arbitrary dielectric 
objects will be provided elsewhere [43|, here we will dis- 
cuss the general theory and present the results obtained 
for objects smaller than the optical wavelength. 

The key term of the total Hamiltonian describing the 
effects of light scattering is i^sc, defined in Eq. (|3T]) . 
First, we are interested in the case in which the cav- 
ity is strongly driven. Then a coherent steady state is 
present in the cavity field, which we explicitly consider 
by decomposing the cavity field into a quantum part plus 
a coherent (classical) part, £^cav(x) £^cav(x) +fcav(x), 
where ^cav(x) is the classical one. A detailed discus- 
sion of this transformation is given in Sec. IV Ai In this 
framework, the scattering term of the Hamiltonian can 
be written as 



free 



Hi 



H^. + Hi 



(34) 



where we have defined: 



HL = -CcCo / C^xEfree(x)fcav(x) 
Jvir) 



^sc = -^c^O / (ix£;free(x)£;cav(x) 
Jvir) 

i^tw = -^c^O / C^x£;free(x)ftw(x). 
Jvir) 



(35) 



Additionally, the displacement of the cavity field also 
modifies the shift term of the Hamiltonian, which 
now reads H^h = -ecCo /(^cav(x) + fcav(x))ftw(x)(ix. 
Whereas the first term is already included in H^^ (since 
in £^free(x) WC integrate over all the electromagnetic 
modes without excluding the cavity mode), the second 
term yields a shift of the trapping frequency as well as of 
the equilibrium position, as discussed in App. IB 31 
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The time evolution of the density matrix describing the 
center of mass motion p is determined by tracing out the 
free modes using a Markovian master equation approach. 
Its utihzation is justified for the fohowing reasons: first, 
due to the fact that the reservoir of free modes of the 
electromagnetic field is very large, the bath density op- 
erators are not significantly changed by the interaction, 
such that one can always assume that its density matri- 
ces remain constant in time pE ^ P£;(0). Second, the 
Markov assumption, stating that the decay of correla- 
tions is much faster than any other time scale in the 
system, Tcorr <^ 'Ts, is fulfilled: the Hamiltonian operator 
contains terms ex J exp(— icc;/ct)(i^/c, with a distribution of 
ujk peaked around cjl, which is the fastest time scale in 
the system. Any correlations in the bath scale with the 
mode frequencies uok and thus decay very quickly. The 
details of derivation will be provided in [43|. In here we 
just report the final result. 



where the first (second) term is rg^(rg^^) and Pt denotes 
the power of the trapping laser. Here, nph is the number 
of steady state photons in the cavity due to the driving 
field. We remark that surprisingly these results are in 
agreement with the ones obtained in the standard the- 
ory of decoherence [5Ql-[52j. Another important remark 
is that the dissipative term V^[p\ does not create any 
mechanical damping of the oscillator, but only diffusion. 
Hence, in this case it might be misleading to discuss the 
mechanical quality factor of the harmonic oscillator. We 
think that it is more appropriate to consider coherent 
times l/Fsc as usually done in the case of ions. The har- 
monic oscillator can oscillate without mechanical damp- 
ing [13 for very long times, however, a quantum state 
prepared in the harmonic oscillator will lose coherence in 
a time scale of l/Fsc- Using typical numbers, see App. [El 
this corresponds to time scales of the order of 0.1 ms for 
nanospheres. 



pit) = i[p, i^oM + ^rn] + >^cav[p] + >^sc [p] + ^m[p]- (36) Finally, H^^ and H^^ also yield an additional force i^m, 



Comparing to the case without scattering, Eq. (|32]) . the 
new terms are the following: first, two dissipation terms 



^sc[p\ = Ai:sc(2apa^ — a^ap — pa^a), 

I>m[p] = r,e [b + b\ [b + b^,p]] . 



(37) 



The first one C^c [p] describes cavity photon losses due to 
events in which cavity photons are scattered out of the 
cavity. This term, which contributes to the decay rate 
of the cavity, is obtained from H^^. For spherical objects 
smaller than the wavelength, a^sc is given by 



(38) 



where we have assumed the sphere to be trapped at the 
maximum slope of the standing wave, kcZ ~ 7r/4, and 
the cavity volume is defined as Vc = 7r/{4:dW^). As a 
check of the theory, one can compare this expression with 
the decay rate one would estimate using the Rayleigh 
cross section aR. With this, the optical finesse is esti- 
mated as J^R = irW^/aR^ and consequently the decay 
rate k.r = caR/{4:Vc). The Rayleigh cross section ne- 
glecting the different polarization of the incoming and 
scattered light (to be consistent with the rest of the ar- 
ticle) is 



CTR 



47r 



(39) 



By plugging ctr into f<iR = caR/{4:Vc) one recovers the 
same expression as derived in the theory, Eq. ([38|) . 

The term V^lp] describes recoil heating due to elastic 
scattering out of the cavity. This is obtained from H^^ 
and H^^. The heating rate Fsc = F^^^ + F*^ has two con- 
tributions from the cavity photons and from the tweezers 
light. For sub- wavelength dielectric spheres, it reads 



; ' 



(40) 



which for subwavelength spheres modifies the trapping 
frequency by cjt ^ Cc;t(l — ec), and the optomechanical 
coupling by go ^ go{l — Cc). This contribution has to 
be added to the optomechanical Hamiltonian i^oM, and 
represents a non-negligible QED renormalization of the 
Hamiltonian due to virtual photon exchange. This QED 
effects will be addressed in Pfianzer et al.^ in preparation^ 
where we will show that higher perturbative terms, e.g., 
corresponding to emission and reabsorption of two pho- 
tons are suppressed by several orders of magnitude for 
small spheres. 



IV. QUANTUM ELASTICITY 

Let us now address the coupling of the center of mass 
(CM) motion mode to other internal vibrational modes. 
One can model the dielectric as an object containing 
constituents, in this case atoms, that are coupled to each 
other by mutual interactions, here modeled by springs. 
The entire nanodielectric inherits N different modes, one 
of them is the center of mass mode; a collective movement 
of all the system's constituents into the same direction. 
The other modes can be described as movements of the 
different constituents relative to each other, mediated by 
the springs. All of these different modes are also coupled 
to each other, which, in turn, infiuences their form and 
lifetime. In principle, one can couple any of these modes 
to light, especially if the object is sufficiently large. We 
are particularly interested in the CM mode in this article. 
We will focus on investigating the infiuence of the relative 
modes, also denoted as vibrational modes, on the center 
of mass mode treating them as a source of decoherence: 
the vibrational modes can in principle take the role of a 
thermal bath and prevent ground state cooling of the CM 
degree of freedom. In order to investigate this source of 
noise, we use an elasticity theory for quantum systems in 
this section. After introducing a field characterizing the 
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FIG. 2: Coordinates used to describe a position 
within an arbitrary dielectric object given by = r + 
02, 03) (u(x) + x), where r denotes the center of mass, 
u(x) a small displacement from the equilibrium position x and 
i?(0i,02,03) the Euler rotation matrix acting on the entire 
object. 



object's deformation, we determine the vibrational eigen- 
modes in Sec. IIV Al Thereafter, we analyze the effect of 
an additional external potential and the induced interac- 
tions between CM and vibrational modes in Sec. IIVBI 
Finally, in Sec. IIV CI we discuss the effect for small ob- 
jects and obtain an effective Hamiltonian by adiabatically 
eliminating the internal modes. 



A. Vibrational eigenmodes 



Let us start by defining the coordinate x', which describes 
a point in the dielectric object. As illustrated in Fig. [21 
this can be written in the most general form as 



x' = r + i?(0i,02,03) (u(x) +x) , 



(41) 



where r denotes the center of mass position. In the 
coordinate system centered at the center of mass posi- 
tion, X is the coordinate describing an equilibrium point 
and u(x) its deformation field. The term i^(0i, 02,03) 
is the Euler rotation matrix with the Euler angles 
01,02,03 that is used to rotate the coordinates x and 
u(x). Note that the center of mass position can be 
defined as r = J d^xp{x.)x.' /{J d^xp{x.))^ with p(x) de- 
noting the system's density distribution. Therefore, 
/ (i^x/)(x)[x u(x)] = 0. In order to guarantee that r 
remains the CM coordinate in case of a vanishing defor- 
mation field, i.e. u(x) = 0, one requires J d?xp{'x)^ = 0, 
and consequently, the deformation field always has to ful- 
fill / d^xp(x)u(x) = 0. 



three-dimensional potential V{x.') reads |53L l5 



jC 



i^(x)x''-nx')-VE(x) 



The elasticity potential is given by 

^e(x) = ^^(Jij(x)ey(x), 



(42) 



(43) 



where 



1 / dui{-x) duj{-K) 



dxj 



dxi 



k 



(44) 



are the elasticity and the stress tensor. The Lame con- 
stants are defined as A = cry[(l + <j)(l — and 
/i = y[2(l + <j)]~^, with (J being the Poisson ratio and Y 
the Young modulus characterizing the elastic properties 
of the material. One can now replace the expression of 
x' in the kinetic part of Lagrangian and obtains 



(45) 



d^x[V{x') + VE{^)], 



where the dots denote time derivatives and li is the ob- 
ject's moment of inertia. We have used that in the kinetic 
part of the Lagrangian, the rotational, vibrational, and 
center of mass degrees of freedom decouple (5^ . 

Let us now determine the unperturbed vibrational eigen- 
modes of the system, that is, the modes obtained without 
considering the potential density V^(x'). In the follow- 
ing, we will assume for simplicity the homogenous case 
p(x) = p, the non-homogeneous case can be incorporated 
easily. Also, we will omit the rotational modes since they 
decouple without the presence of the external potential. 
Let us first derive the Hamiltonian by defining the CM 
momentum as pi = dC/dvi and the momentum density 
as Vi{yi) = dC/dui{yi)^ leading to 



H. = 



P_ 

2M 



[v(x)]^ 
2p 



V^e(x) 



(46) 



One can determine the vibrational eigenmodes by sepa- 
rating variables in the corresponding equation of motion 
for u(x, t), which reads [53| 



pu(x, t) = /iV^u(x, t) 



-AV[V-u(x,0]. 



(47) 



Here, u(x, t) can be separated into transversal and lon- 
gitudinal oscillation modes, u(x, t) = u^(x, t) + U|| (x, t), 
where V • u^(x, t) = and V x U||(x, t) = 0, and either 



The system's Lagrangian in the presence of a general open or periodic boundary conditions can be used. The 
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longitudinal modes describe compression waves propa- 
gating at velocity c\\ = [(A+2/i)/p]^/^ and the transversal 
modes torsion wave propagating at c± = [/i/p]^/^. In the 
following we will only consider the longitudinal modes 
along the cavity axis. By expanding the elasticity field 
along the cavity axis for these eigenmodes u^{z) (which 
are normalized as Jy d^xu^{z)u^{z) = (^nm"^)? one has 



u{z,t)=^u'i{z)Qn{t) 

n 

V{z,t) = Y,ul{z)Pn{t), 



(48) 



where Pn{t) = pQnif)- By plugging this decomposition 
into the Hamiltonian Eq. ([46|) , one obtains after some 
algebra 



(49) 



where the frequency of the internal modes is given by 

n 2 



2 _ A(l-a) 



[ 

Jv 



d 

dz 



(50) 



The eigenmodes u^{z) have to be chosen accordingly to 
the geometry of the object. We will discuss the specific 
form of the mode and the value of the parameters in 
Sec.lED 

At this position, it is straightforward to perform a canon- 
ical quantization of the eigenmodes Qn, by considering 
them as operators fulfilling the canonical commutation 
rules [Qn^Pm] = i- As already done in the previous sec- 
tions, the momentum operator of the CM will also be 
quantized with the external harmonic trap. 



B. Effect of the external potential 



The external potential V{x.') can in principle effect a 
coupling between the rotational, the center of mass, and 
the vibrational degrees of freedom. In case of a purely 
isotropic harmonic potential, it can be easily verified that 
the coupling vanishes. On the other hand, for arbitrarily 
shaped objects, the external anharmonic part of the po- 
tential effects some coupling between all degrees of free- 
dom. In the following we assume spherical objects, for 
which the direct coupling between the CM and the rota- 
tional degrees of freedom vanishes. Even in the case of a 
prolate spheroid, the coupling is negligible [2]. For spher- 
ical objects, there is only an indirect coupling between 
the CM and the rotations, mediated by the vibrational 
modes, which is negligible and will be omitted hereafter. 
Therefore, with these assumptions one can consider the 
center of mass mode to be decoupled from the rotational 
motion, and we consequently omit the rotational modes 



in the rest of the section. One can then focus on the 
one-dimensional case derived in the previous section by 
only considering the longitudinal modes. 

The total Hamiltonian, including the external potential, 
is hence given by 



H = H, 



+ / d^xV{z') 
Jv 



(51) 



Assuming that the deformations u{z) are small and that 
the object is trapped at r ^ 0, one can expand V{z^ = 
z + u{z) + r) to second order in r and u{z)^ which leads 
to 

H=Ho^r [ d^xV'{z)^— [ d^xV"{z) 

Jv 2 Jy 

+ d^xul{z)ul,{z)V"{z) 

(52) 

+ y,Qn d^XU^^{z)V\z) 

^rJ^Qn f d^xu'^{z)r\z), 

where the primes denote spatial derivatives. By recalling 
that the external potential is, in our case, given by the 
light-matter interaction term Eq. (|24|) . that is V^(x') = 
— eceo[£^(x')]^/2, one can understand the terms appearing 
in Eq. ([52]) as follows: 

1. The term r jy d^xV {z) yields the optomechanical 
coupling of the center of mass mode as described 
in App. [B2l 

2. The second term r'^ jy d^xV" {z) /2 describes the 
harmonic trap of CM mass given by the optical 
tweezers, as described in App. IB II 

3. The term QnQm jy d^xul{z)ul^{z)V" [z) /2 gives 
a correction to the harmonic trap for the inter- 
nal modes as well as a coupling between internal 
modes. 

4. The first new interesting term is 
Qn jy d^xu^^{z)V' {z)^ which describes an op- 
tomechanical coupling between the internal modes 
and the cavity field. 

5. Finally, the most relevant term for our purposes 
is rQn jy d^xu^{z)V" {z)^ which describes the cou- 
pling between the vibrational degrees of freedom 
Qn and the center of mass mode r. 

Taking into consideration these terms, one can now write 
the center of mass mode as r = ^o(^^ + where zq is 
the ground state size, as used in Sec. lIIH and the internal 
modes as Qn = qo,n (cn + 4)^ with (7o,n = (2Mcj^)"^/^. 
Note that, due to the additional external tr aping with 
frequency cjt, the effective vibrational frequencies are 
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changed to uj'^ = {uj^ + ^n)^^^ ^iii omit the prime 
hereafter). The new part that has to be added to the 
total Hamiltonian i^tot, see Eq. (|29]) , which takes into 
account the presence of internal modes, is given by 

He = ^ujncicn ^^gn{a,a'^){cn + 4) 

n n 

CO 

+ + 4)(Cm + cl^) (53) 

n,m 

oo 

+ Y.^r^{Cn + ci){b + b^). 

n 

The coupling between the cavity field (which de- 
pends on the cavity mode a) and the vibrational 
modes is given by gn{a,a^) = qo,n fy d^xV\z)u^{z). 
The coupHng between the internal modes is = 
qo,nqo,m fYd^xu^{z)ul^{z)V{z)/2. Finally, the cou- 
pling between the CM mode and the vibrational modes 
is given by 

7n = zoqo,n I d'xV"{z)ul{z). (54) 

JV 

Summing up this subsection, we have derived the quan- 
tized Hamiltonian describing the coupling between the 
CM and the vibrational modes in the presence of an ex- 
ternal potential density. It can be shown that for a har- 
monic external potential, the CM mode is decoupled from 
the internal ones since V"{z) is constant and by recalling 
that jy d^xu^{z) = 0, one obtains 7^ = 0. In the next 
section, we estimate the order of magnitude of the pa- 
rameters for objects smaller than the optical wavelength 
in the presence of the anharmonic potential given by the 
standing wave. 



C. Sub- wavelength spheres 

First of all, let us estimate the order of magnitude of 
the internal vibrational frequencies, see Eq. ([5Q|) . for the 
case of a sphere of radius R. To get an estimation of 
the order of magnitude, for simplicity one can just use 
the eigenmode u^{z) = cos(knz) with kn = n7r/(2i?), ob- 
tained for a cube of length 2R and with open free periodic 
boundary conditions. Then, using typical values of the 
Young's elasticity module Y and the Poisson constant a 
(see App. [E|) . the vibrational frequences are of the order 
ujn ^ 10^^ Hz {uon ^ nc\\/R). Note that comparing this to 
the typical values of the CM frequency ujt ~ 10^ Hz, the 
internal frequencies are five orders of magnitude larger 
for objects of the order of 100 nm. 

This large difference in frequencies between the CM 
modes and the internal modes enable us to adiabatically 
eliminate the vibrational energy levels. It can be shown 
that this approximation is justified by solving the equa- 
tion of motion for the CM and vibrational operators by 



applying Laplace transformations. The solution obtained 
in this way contains parts oscillating at frequencies 00^ 
and cjn, where all terms oscillating at uOn are suppressed 
by a factor uot/uon ^ 1- Thus, it is well-justified to neglect 
these terms and to perform an adiabatic elimination. One 
can perform this by eliminating the vibrational levels on 
top of the steady state, yielding the result that the only 
effect is a shift of the trapping frequency of the CM mode 
given by 

where {c\cn) is the occupation number of phonons in 
the vibrational mode n. By plugging in typical num- 
bers, one gets a correction to the trapping frequency of 
{uo[ — uot)/oot ^ 10~^^, which shows that the CM mode 
is decoupled from the internal modes for objects smaller 
than the optical wavelength. 



V. LIGHT-MECHANICS INTERFACE 

One of the most fascinating perspectives of quantum op- 
tomechanics is the possibility to prepare superposition 
states of objects containing billions of atoms, and there- 
fore, to test quantum mechanics at larger scales. Already 
in the early days of this research area, several groups 
proposed to create non-classical states of a movable mir- 
ror (55l-[57|. The idea behind these proposals is to use 
the optomechanical interaction to entangle a small quan- 
tum system with the macroscopic object. By observing 
the state of the small quantum system, the creation and 
loss of the non-classical state in the macroscopic system 
can be monitored. This idea was also used in [58], where 
the coupling between a micromechanical resonator and a 
Copper box is proposed in order to prepare entanglement 
between the quantum system (Copper box) and the can- 
tilever. We remark that an experiment has been recently 
reported in [59], where coherent control of single phonon 
has been achieved in a high-frequency micromechanical 
oscillator. In Marshall et al. [21] (see also [22|) a scheme 
to prepare a superposition state of two distinct locations 
of a mirror through the optomechanical interaction with 
a single photon has been proposed. All these ideas pose 
a major challenge to an experimental realization mainly 
due to the following reasons: (i) the coupling between the 
small quantum system and the macroscopic mechanical 
system is not strong enough and (ii) the mechanical sys- 
tem suffers from its fast decoherence due to the thermal 
contact. 

In this section, we show a possible way to circumvent 
these two restrictions. We propose two protocols to 
strongly couple a no n- Gaussian light state to a mechani- 
cal object. This is achieved by using a driving field which 
enhances the interaction into the strong-coupling regime 
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(the interaction time has to be faster than the decoher- 
ence times). This enhancement of the optomechanical 
couphng by the driving field was suggested in [26|, li^ 
and experimentally observed in Then, on top of 

the driving field, which is red-detuned, a quantum light 
state is sent into the cavity which is transferred to the 
mechanical system by the strong coupling. This idea has 
been introduced in [l| (see also [61, 62]). Additionally, 
we propose an alternative protocol that uses the weak 
coupling regime to prepare non- Gaussian states. These 
protocols, which can be applied to general optomechan- 
ical systems, are ideally suitable for optically levitating 
nanodielectrics, since they do not have a thermal con- 
tact [H, [11, and thus possess longer coherent times. 



These general light-mechanics interface protocols allow 
us to prepare non-Gaussian states by using a Gaussian 
Hamiltonian. Their key ingredient is that one uses non- 
Gaussian input states (similar ideas have been used in 
the context of quantum computation [63|, Q). Hence, 
these protocols represent an effective and simple way to 
produce non-linearities in optomechanical systems, a goal 
that is intensively pursued (see for instance (Zlj). 

Finally, we remark that in case of a levitating object light 
scattering yields decoherence of the mechanical state 
with a rate given by Fsc- For sufficiently small objects, 
this can be made much smaller than Hi. In the following, 
where we are interested in designing the protocols, we will 
neglect the effects of light scattering (more precisely, the 
term Hsc) by assuming that the protocols can be realized 
on a time scale much shorter than l/Fgc- For other op- 
tomechanical setups, decoherence in the mechanical sys- 
tem could be incorporated easily into the protocols. 



This section is organized as follows: first, in Sec. IV Al we 
transform the total Hamiltonian of the system in order to 
account for the driving field of the laser. We then present 
three different protocols in Sees. IV B] IV C[ and IV Dl 



A. Driving field: displaced frame & initial state 



In this section we transform the Hamiltonian J^tot, see 
Eq. (|29|) . into Hl^^ in order to incorporate the driving 
field of the cavity and keep a close structure of the Hamil- 
tonian. This allows us to describe the quantum states on 
top of the steady state which will be used in the protocols. 
Throughout the article we will use either the original 
frame, in which states are described according to i^tot, 
or the transformed or displaced frame, in which states 
are related to i^tot- Then, we describe in both frames 
the form of the total initial state that one obtains after 
coohng the mechanical oscillator to the ground state. 



1. Displaced frame 



In this section we will perform the standard transforma- 
tion [26, 28, 44] done in quantum optomechanics in order 
to shift the coherent part of the states obtained when 
driving the cavity with a laser. However, in contrary to 
what is usually done, here we also need to displace the 
output modes since we use them in the light-mechanics 
interface. 

First, one moves the cavity and the output field to the 
frame rotating with the laser frequency cul- This is de- 
scribed by the unitary operator 



Ur{t) = exp 



-IUJl 



aQ{uj)ao{uj)duj 



(56) 

To ease the notation, after this transformation we rede- 
fine the ao(co') and j{oo) such that ao(cj) = ao(a; + cjl), 
and 7(cj) = ^{uj -\- ujl)- The total Hamiltonian (ignor- 
ing the scattering part Hsc and the shift H^h which is 
discussed later) reads 

/CO 
ujal{uj)ao{uj) 



J — c 



^{uj){a^ao{uj) - H.c), 



(57) 



where Aq = uJc — ojl- Then, one displaces the cavity field 
with the displacement operator Da{cy)^ the mechanical 
field with D}){/3), and the output modes with Dout{a^), 
that is, 



Dl{a)aDa{a) = a 
Dl{f3)bD,{/3) = b^ 



(58) 



After applying this transformation to the Hamiltonian, 
one fixes a, /3, and a^^, such that the terms in the Hamil- 
tonian that have only one creation or annihilation opera- 
tor vanish. This corresponds to solving the following set 
of equations: 



Aoa + 2goaf3 + i 
uJtf3^go\a\^ =0, 
ujal{uj)auj - 



j{uj)al{uj)a = 0, 



(59) 



which have the solutions 



lA 



9om 

OJt 

( 

V7(0) 



(60) 



7ra7(0) I 5{uj) + ia'^{uj)V{uj 
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Here, A = Aq + 2^0/5, and ao = ^1^/7(0), where 
Ql = yJ'^Pc^l^L^ Pc being the laser power. The sym- 
bol V denotes the principal part, and we have used that 
7^(cj) ^ tz/n in a finite region around uj = [42| in or- 
der to perform the integral V J^^(jO~'^duj = 0. In the 
next subsection Sec. IV A 2[ we show how to obtain the 
expression of from a more physical perspective. 

To sum up, the transformation applied to the Hamilto- 
nian can be defined as P = Dout{aoj)Di){l3)Da{(^)^ and 
the transformed Hamiltonian is given by 



LC5 



where 



(61) 



(62) 



is the enhanced optomechanical Hamiltonian, and i^LC 
is transformed into 



/oo 
uja\{ijj)aQ{(jj)dijj 
-UJl 

/oo 
'y{uj){a^ ao{uj) — I{.c.)d(jj. 



(63) 



Note that Eq. ()6T|) has the same structure as Eq. ()57j) 
with the only replacement Aq A, and goa^a{b^ -\-b) ^ 
g{a^ + a){b^ + b). We have defined g = ^o|<^|, and 
^ = arg(a), and we have redefined the a {ao{uj)) opera- 
tors as = ae~^^ i^oi^) = <^o('^)e~^^) (we omit the tilde 
hereafter). A crucial remark is that the optomechanical 
coupling g is enhanced by a, which is the square root 
of the mean number of photons inside the cavity in the 
steady state. This will allow us to reach the strong cou- 
phng g ^ n (where n is the decay rate of the cavity) in 
the light-mechanics interface. 

We remark that in case of using levitating objects, the 
shift to the trapping frequency as well as the shift in 
the equilibrium position, given by the Hamiltonian i^sh, 
should be taken into account in the i^oM Hamiltonian. 
As discussed in App. IB 3[ this would imply to change 
the trapping frequency to uot ^ oot ^ and the dis- 
placement of the cavity mode to {3 /3 + ^sh/^, where 
cjsh and ^sh are given in App. IB 31 However, to keep the 
section in a general form, so that it can also be applied 
to other optomechanical systems, we will omit this effect 
hereafter. 

The transformed Hamiltonian can now be written in the 
interaction picture (assuming that the free part is Hq = 
utb^ b -\- Aa^ a -\- ujal{uj)ao{uj)duj) as 



/oo 
7(w)(atao(u;)e'('^*-'^*) - H.c.)dw. 



(64) 



Now, by choosing a red-detuned driving A = cj^, one 
can perform the rotating wave approximation (valid at 



^ g)^ and obtain the beam-splitter interaction form 
of the total transformed Hamiltonian in the Schrodinger 
picture 

Hl^^ = ujtia'^a + 6^6) + gia'^b + a6^) + Hie (65) 

Analogously, one can consider a blue-detuned driving 
A = —ujt in order to get the two mode squeezing in- 
teraction Hamiltonian: 

H^^^ = - ujtia'^a - b^b) + g{o)b^ + ab) + Hi^ (66) 

These two types of interaction will be used in Sec. |Vl 
to design different protocols in the light-mechanics inter- 
face. 



2. Initial state 



All the protocols that we shall discuss in the next section 
assume that the initial state is the ground state cooled 
by the red-detuned field (A = cjt). As discussed in the 
previous section and in App. [Dl this state is given by 

/CO 
I>outK)dw|0)out = I>|00n), 
-UJl 

(67) 

where "b (a)" labels the subspace of the mechanical mode 
(cavity mode), "out" the subspace of the output modes, 
and the vacuum state for the output modes. The dis- 
placements a, /3, and are defined in Eqs. ([6Q|) . 

Note that |in) is an eigenstate of the total Hamiltonian 
i^tot, see Eq. (f57j) . This can be trivially demonstrated 
by using that V^HtotV = H^^^ (for the red-detuned case 
Eq. ([65]) ). and that Hl^^\00^) = 0, since then one has 

i^totlin) = VV^HtotV\00n) = VH^JOOQ) = 0. (68) 

The state |in) (reading \00Q) in the displaced frame) will 
be considered as the initial state upon which the pro- 
tocols are designed using either the beam splitter inter- 
action Eq. ([65]) or the two mode sequeezing interaction 
Eq. ([66]). 



B. Reflected One-photon 



In this section, we will present a protocol which strongly 
couples a one photon state to the mechanical motion of 
the oscillator. This protocol is general and can be applied 
to various optomechanical systems. Let us remark that 
it has already been introduced by some of the authors in 
[H and that related ideas have been reported in [6l|, |62| . 
In this section we will provide a thorough analysis. In 
particular, we develop a formalism to solve the input- 
output formalism in the Schrodinger picture in order to 
be able to describe the final state of the protocol. 
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Let us start by sketching the different steps of the proto- 
col: 

1. Cool the mechanical motion to the ground state by 
the red-detuned driving field. 

2. Keep the strong driving field switched on such that 
the beam-splitter interaction is induced inside the 
cavity. 

3. Impinge the cavity with a resonant single-photon 
state, sent on top of the driving field as a result of 
parametric down conversion followed by a detection 
of a single photon j45| . 

4. When impinging the cavity, part of the field is re- 
flected and part transmitted (65| . 

5. The beam-splitter interaction Eq. (|65]) caused by 
the red-detuned laser, swaps the state of light in- 
side the cavity to the state of the mechanical mo- 
tion. 

6. By tuning the width of the light pulse appropri- 
ately, one finds that at time th, one has a maxi- 
mum mean number of phonons of 1/2 in the me- 
chanical system. At that time, the driving field is 
switched off. Then, the entangled state |£^)out,b ^ 
|0)out|l)b + e'^|i)out|0)b is prepared. Here out(b) 
stands for the reflected cavity field (mechanical mo- 
tion) of the system, and |0(l))out is a displaced vac- 
uum (one photon) light state in the output mode 
of the cavity Aout- The phase ^, given by the light- 
mechanics interaction, is always fixed. 

7. At a later time, once the reflected photon is far 
away from the cavity, a balanced homodyne mea- 
surement of the output mode is performed. The 
motional state collapses into the superposition 
state |^)b = co|0)b + cie^'^|l)b, where the coeffi- 
cients Co(i) depend on the measurement results. 

In the following we will analyze carefully the important 
steps of the protocol. In the shifted frame, the initial 
state (according to Sec. IV A 2]) . consisting of a photon on 
top of the ground state of the mechanical oscillator, is 
given by 



/oo 
<l>Uu)al{co)\00n), 
-UJL 



(69) 



where shape of the photon pulse which is 

assumed to be Gaussian 



0in(^) 



1/4 



(70) 



Here, x^^ is the position from which the pulse has been 
sent (it is considered to be large, Xin ^0). A = 
— ool = oJt is the detuning, which shows that in the 



non-rotating frame the pulse is centered at the resonance 
frequency of the cavity. Note also that one can express 
the mode function in position space by the Fourier trans- 
form (^in(x) = / (icj0in(cj)e^^^/V27r. 

The time evolved state with the beam-splitter interaction 
Eq. (|65]) . I'^(t)) = exp[—ii^[Q^t] 1^/^(0)) can be expanded in 
the following basis. 



\^(t)) =cb(t)\im) ^ ca(t)\{)m) 



f 

J — c 



t)al{uj)duj\mQ.). 



(71) 



The time-dependence of the coefficients can be obtained 
using the Wigner-Weisskopf formalism. By using the 
Schrodinger equation, one obtains 

Ch{t) = -\ujtci,{t) - \gca{t), 

/oo 

c{uj^t) = —iujc{uj^t) — 'y{uj)ca{t). 

(72) 

This system can be further simplified by formally solving 
the differential equation corresponding to c(a;,t), plug- 
ging it into the equation for Ca(^), and by using the ap- 
proximation j{uj) ^ 7(0) = \f^J^. One gets (analagous 
manipulations have been explicitly done in App. [D|) : 



Chit) 
Ca{t) 



c{uj^ t) 



-\ujtCh(t) - \gca(t), 
-(\ujt + i^)ca{t) - igcb{t) 

L 



(73) 



This system of differential equations can be solved by 
using that Ca(0) = 05(0) = and c(cc;,0) = (j)^^{uo). After 
some effort, one obtains 



Ca(t) = fp_{t - T)4>UT)e-'^^^'--U' 

Jo 

Cb{t) 

The functions q{t) and p-{t) are defined by 



t 





(74) 



P±{t) 



-Ktl2 



cosh(xt) ± — sinh(xt) 



-i^e-"*/2 sinh(xt), 

X 



(75) 



where x — \/ ^^/4 — g'^. One can now plot the mean 
number of phonons ni){t) = |c5(t)p and photons na{t) = 
|ca(t)P, see Fig. [3] for some parameters given in its cap- 
tion. Note that at t = th^ where 



arccos(Ai:/2^) 



(76) 
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Fig. 3: Input-output dynamics after sending a one-photon 
pulse centered at xin = 5//^ from the cavity at t = 0. A 
Gaussian pulse of width a = 5.6/^ is used. We plot the mean 
number of phonons in the mechanical system nh(t) — \ch{t)\^ 
(red solid line) and the mean number of cavity photons 
na{t) — \ca{t)\^ (blue dashed line). We consider the strong 
coupling regime g — and tune the width of the pulse so 
that the maximum mean number of phonons is 1/2 (dotted 
grey line) at t = t^. 

the mean number of phonons 725 is maximal. By tuning 
the width of the initial pulse, one obtains that Ca{th) ^ 
and \ci){th)\ ^ 1/a/2. In this case, the total state at th is 
given by 

/oo 
c{uj,th)al{uj)duj\OOn). 
-ojl 

(77) 

This is an entangled state between the ouptut photon 
mode, described by the pulse shape c(a;, th)^ and the me- 
chanical phonon mode. In the non-displaced frame, the 
state at th is described by \ilj\th)) = V\il){th)) - 

At t = th the driving field is switched off. However, at 
this time, there is still a large number of photons |ap 
present inside the cavity. They will leak out of the cav- 
ity reducing the classical force that they were exerting 
on the mechanical system, which is described by the dis- 
placement of the mechanical system, p. In order to com- 
pensate this effect, one could move the center of the trap 
muj^{x — Xt{t))'^/2 accordingly, which yields a force term 
— mco'^Xt(t)xo(6 + 6^), in order to keep the ground state of 
the harmonic oscillator. Another effect of this leaking out 
of photons is that the coefficient ci){th) will be decreased 
at some later time. Note however, that one could send a 
pulse that generates \c}){th)\ > 1/a/2 such that, after the 
decrease due to the leaking out of the coherent photons, 
one obtains \ci){t > th)\ = l/\/2. The discussion on how 
to compute and estimate this effect is done in App. ID 31 

Here, we just approximate the state dX t^ th as 

m)) = C6(t/.)e-^"*(^-^^)l^out|101^) +I^out4utt|0017), 

(78) 



where I^out only displaces the output modes with a^. 
Also, we have defined the output mode of the cavity as 

/oo 
(l)out{oJ,t)al{uj)duj, (79) 

where ^out(^,^) = c{uj,th)e~'^^^~'^''\ Note that the dis- 
placement is only in the output modes since the photons 
inside the cavity, and the consequent radiation force into 
the mechanical object, are not present at times t ^ th 
since the driving field is switched off. 



1. Measurement of the output mode 

The final step of the protocol is the measurement of the 
quadrature of the output mode Aout,t, that is 

-^out,t = ^lut,t + ^out,t- (80) 

This measurement consists in integrating the signal of a 
continuous measurement with the mode shape given by 

0out(^,^)- 

More generally, the output operator Aout = 
S'^ojL can be written as a combina- 

tion of mode operators at position x by using 
ao(co') = dxe~'^^^ao{x)dx/\^^^ leading to 

Aont = / 4>{x)aQ{x)dx. (81) 
Jo 

Note that now the mode ao{x) can be measured at the 
position X = Xd of the detector at time t, by the relation 
doixd^t) = ao{x = Xd — t^O). Then, by a continuous mea- 
surement of ao{xd^ t), one has access to the measurement 
of all ao{x) and consequently, also to Aout by integrating 
the signal over (j){x) (note that Aout is a linear combina- 
tion of the independent modes ao{x)). 

After the continuous measurement, let us assume one 
obtains the value Xout- Then, the superposition state in 
the mechanical object, given by 

|*> = -i=(co|0)6 + ci|l>fc), (82) 

is prepared, where Co(i) = (^out|l(0))- The measurement 
of the quadrature poses an experimental challenge. If we 
define the two orthogonal states |±)o = 1^^) ± aI^^^\Q) 
and their displaced states |±) = I)out|i)o7 one obtains 
that the mean value and fluctuations of Xout are given 

by 

(-^out,t)± = Oix -\- (Xout,t)±,07 

/2\ 2 / \ /2\ \^oO J 

\-^out,t)± = ^x^ 2aa,(Xout,t)±,o + \-^out,t)±,0, 

where we have defined = L>Jut,t-^out,t^out - ^out,t- 
Thus, (AXout,t)± = (AXout,t)±,o- This shows, that from 
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the theoretical point of view, the two displaced states 
|±) are as distinguishable as the non displaced ones |±)o- 
From the experimental point of view, the problem is that 
the signal to noise ratio in a balanced homodyne mea- 
surement is too low. Although the displacement ol^ can 
be computed by using 0out(^,^) and (see Eq. ([60]) ). 
the final expression is not very illustrative. Instead, in 
App. ID 2\ we will analyze the problem of the measure- 
ment of the output field when a photon on top of the 
coherent field was prepared inside the cavity. This gen- 
eral problem is more enlightening and the conclusions 
apply directly to the reflected-photon protocol, that is, 
that the displacement of the output field is of the order 
of a - 10"^. 

In order to circumvent this experimental challenge we en- 
visage the following ways out: i) subtract the coherent 
part by destructively interfering a coherent beam with 
the same phase; ii) use an optomechanical system where 
the detuning between the resonant photon and the red- 
detuned driving is much larger (since A = cjt, this would 
correspond to a large- frequency mechanical oscillator). 
This must be done without compromising the strong cou- 
pling requirement which is based on the enhanced cou- 
pling g = I^^lI^o/VA^ + iii) use a scheme similar to 
the one proposed in [62|, where the photon is sent in 
the dark port of an interferometer; iv) design a scheme 
where the light pulse is perfectly absorbed in the cavity 
and therefore no measurement is needed. In the following 
section we present a protocol which follows solution iv). 



C. Perfect absorption 



Here we present a protocol which circumvents the chal- 
lenging step of measuring the displaced output mode in 
the reflected one-photon protocol, as discussed in the pre- 
vious section (this protocol was announced in [1]). The 
goal is to perfectly absorb the light pulse, which is in a 
non- Gaussian state, into the cavity, and therefore trans- 
fer it to the mechanical system. This is achieved by using 
a time modulation of the optomechanical coupling g{t)^ 
which can be implemented by varying the intensity of the 
driving field. Similar ideas have been proposed in the 
context of quantum communication |66| , and recently in 
quantum optomechanical transducers [13 • 



1. Time- dependent displacement 



In this section, one cannot use the beam-splitter inter- 
action (|65]) . since the laser intensity is time-dependent. 
Therefore, a time-dependent displacement has to be per- 
formed carefully. Let us start with the basic Hamiltonian 
in the non-displaced frame Eq. ([57]) . From there, one can 



derive the evolution equations for a, 6, and ao{uj) 

/oo 
duo'^{uj)aQ (cj, t). 



b — —iujtb — igoa^a^ 

do(co', t) = — icjao(cj, t) — j{uj)a. 



(84) 



The tool that will be used in this section is a time- 
dependent driving field at the laser frequency a; = (in 
the rotating frame). This can be incorporated by apply- 
ing the following displacement to the output modes 



ao{uj,t) ao{uj,t) - a -QL{t)6{uj). (85) 
V 

By formally integrating the equation of do(cj,t), and us- 
ing the Markov approximation ^{cj) ^ y^^^T^^ ^^le system 
reads 



d = -(iAo + tz)a — igoa{b^ + 6) + 1]l(^) + \/2^ain(t), 
b = —iuJtb — igoa^a^ 



do(cc;, t) = —iujao{uj^ t) — '^{uj)a + \ j —QL{t)^{^)^ 



(86) 



where we have defined the so called input operator as 



ain(t) = (27r)-i/2/da;ao(a;,0)e 



Next, we perform 



the following time dependent displacement 
a{t) ait) +a(t), 

b{t)^b{t)^m. 



(87) 



and choose a(t) and /3(t), such that the constant terms 
in the equations for a{t) and b{t) vanish, that is 

a = -(iAo + A>:)a - i^o<^(/5 + Z^'') + ^L, 
/3 = -\ujtl3 - \gQ\a\^. 

Then, we perform the following changes of variables 



a{t) ^ a{t)e-'^^\ 
b{t) 6(t)e-^^*^ 

90 



(89) 



{g{t) is real) and perform the RWA considering the red- 
detuned case Ao = uJt- Putting all these things together, 
Eqs. ([86]) read 



d = -K.a - \g{t)e'^b + V^ain(t)e^^o^ 
b = —\g{t)e~^^a^ 

do(cJ, t) = -icjao(cj, t) - 7(cj) [ae"^^°* + a{t)\ (90) 
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We have neglected the smah terms (not proportional to 
a) —igoa{b + b^) and —igoa^a in the equation of motion 
for b. We have also neglected the term —igoa{l3 + P"^) 
in the equation for a. This term, which is smaller than 
— i^foaAo, makes the equation decribing the shape of g{t) 
(to be derived below) much more difficult to solve and is 
therefore also neglect since it does not change the physics 
of the problem. 



Finally, note that equations (f88|) give the solution of 
the time-dependent laser amplitude ftL{t) such that the 
time-dependent coupling g{t) is implemented. In the next 
sections, we derive the pulse g{t) for which any light 
state is absorbed into the cavity and therefore perfectly 
mapped into the mechanical system. 



2. Condition for perfect absorption 

The formal condition for perfect absorption can be de- 
rived as follows. After the transformations are made, the 
evolution equation for ao{uj^t) reads 

do{uj^t) = — iujao{uj^t) — ^{uj)ae~^^^^ 

(91) 

y Ki 

By formally integrating this equation for the initial con- 
dition t = 0, as well as for the final condition t = ti^ and 
subtracting these two solutions after integratin g ove r cj, 
one obtains (using the approximation 7(0;) ^ \f^J^ and 
that 1]l(0) = VtLitx) = 0) 







1 r 

- V2k {a{t)e 



iAo* 



a{t)). 



(92) 



This is the so called input-output relation [42|, which 
relates the output field (the second term containing the 
ao{uj^ti) modes) with the input field ain(t), the quantum 
field from the cavity a(t), and its coherent part a{t). The 
condition for perfect absorption is that the output field 
only contains the coherent part from the cavity, that is 



1 



e-^^it-t^){ao{uj,ti)) = -V2i^a{t). (93) 



- 27r 

With this condition, Eq. ([92]) reads 



(94) 



One can now plug this condition into the Eqs. (j9Q|) and 
obtains 

{d(t))=K{ait))-igit){bit))e'i, 

m) = -ig{t){a{t))e-'i, 
which can be further simplified to 

7l{t)g{t)-fi{t)9{t)+9^{t){a{t)) =0, 
where ri{t) = K{a{t)) — {a.{t)). 



(95) 



(96) 




Fig. 4: Perfect state transfer of a |0) + |1) photonic state by 
sending a Gaussian light pulse of width a = 2k/ 3 from a 
distance Xin = IOk (c = 1). We plot the time modulation of 
g{t)/K, (solid red line) and {b{t)) (dashed blue line). After the 
modulation, when g{t) = 0, one obtains that (b) = 1/2. This 
shows that the superposition state has been mapped to the 
mechanical system without requiring a measurement. 



3. State-independent pulse 



The solution of the equation Eq. ([96]) yields the optome- 
chanical pulse g{t) necessary to perfectly transmit a light 
state into the mechanical system. In order to obtain a 
state-independent solution, we will assume that a coher- 
ent state with phase is sent to the cavity and will show 
that the solution does not depend on ag. Therefore any 
linear combination of coherent states (and therefore any 
state since they form a complete basis) will be perfectly 
transmitted to the cavity with the pulse g{t). 

The initial state is assumed to be 



|^(0)> = exp 



(j)^^{uj)al{uj)duj — H.c. 



(97) 

where shape of the pulse. One can then 

obtain that (ain(t)) = ^^(^^^(t), and using Eq. ([94]) . 
{a{t)) = as4>tn{t)e'^''^/V2i^. Then, Eq. reads 

[/./i(t)-/i(t)]^(t)-[/./i(t)-/i(t)]^(t)+/i(t)^^(t) =0, (98) 

where fi{t) = ^J'^(t)e^^°^ This is the main result of the 
section since its solution yields the time-dependent cou- 
pling g{t) for perfect mapping of any light state into the 
mechanical system, since it does not dependent on the 
coherent phase a^. In Fig. [H the solution g{t) is plotted 
considering (/)in(^) to be the same Gaussian pulse used in 
the reflected one-photon protocol, see Eq. ([7Q|) . 

As an example, let us assume that one wants to transfer 
a photon in a superposition state described by 



1 1 



(99) 
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FIG. 5: Schematic representation of a light- mechanics inter- 
face of teleportation in the bad cavity hmit. The cavity is 
driven by a blue-detuned laser which induces a two mode 
squeezing interaction between the cavity mode and the me- 
chanical mode. Being in the bad cavity limit k > g, the cavity 
photons, which are in a two mode squeezed state |TMS) with 
the mechanical phonons, rapidly leak out of the cavity. The 
output field is combined in a beam-splitter together with the 
non- Gaussian state to be teleported \ip)e- A measurement of 
the output quadratures would realize the Bell measurement 
required for teleportation [gJ- 



In Fig. [Hthe mean value of b{t) is plotted using the g{t) 
solution obtained for the Gaussian case. As expected, 
{b{t)) attains the value 1/2, showing that the superpos- 
tion state (|0) + |l))/\/2 has been prepared. 

To sum up, this protocol enables one to perfectly map 
any state of light into the mechanical system without 
performing any measurement, merely by using a smooth 
modulation of the optomechanical coupling. 



D. Teleportation in the bad cavity limit 



The two previous protocols required the strong coupling 
g ^ n. In this section we provide a protocol, announced 
in [H, to map non- Gaussian states which is also appli- 
cable in the bad cavity limit n > g. The key ingredient 
of the protocol is to drive the cavity with a blue-detuned 
field in order to obtain a two mode squeezing interaction, 
see Eq. ([66|) . The two mode squeezed state is then pre- 
pared by the optomechanical coupling between the me- 
chanical mode and the cavity mode, which rapidly leaks 
out of the cavity. The output mode of the cavity, which is 
in a two mode squeezed state with the mechanical system, 
can be used as an entanglement channel to teleport [g^ 
a non-Gaussian state of light from outside the cavity into 
the mechanical system (see Fig.[5]for an illustration of the 
protocol). This protocol has first been introduced as an 
interface between quantum dots in optical cavities (68| . 
In reference [68|, a detailed discussion of the protocol is 
provided, which applies to our optomechanical setup in 
complete analogy. Thus, we will only summarize and 
remark the important aspects of the protocol here. 

Using the Hamiltonian Eq. ([66]) , one can obtain the equa- 



tions of evolution for a and 6, 

a{t) = -(iAo + n)a{t) - igh\t) + V^ain(t) 
hit) = -luothit) -\ga\t). 



(100) 



One can now move to the interaction picture {a{t) = 
ai{t)e~^^^'^ and h{t) = 6/(t)e~^'^*^), and by considering 
the bad cavity limit ^ g) one can adiabatically elim- 
inate a/(t), by setting d/(t) = 0. One obtains 



hi{t) = !^bi{t) 



^9 



lenity 



■iAot 



(101) 



Formally integrating this equation and using the initial 
conditions (a-^(t)ain(t)) = {b\{0)bi{0)) = (the mechan- 
ical initial state is assumed to be in the ground state) 
yields 



{b\{t)bj{t)) 



1. 



(102) 



This can be used to obtain the squeezing parameter r 
of the entangled state, which will provide the fidelity of 
the teleportation scheme. As proved in [68], the output 
mode of the cavity and the mechanical system are in the 
two mode squeezed state |TMS)5,out, defined by (in the 
displaced frame) 



|TMS)b,out = S (re'"^) |00) 

= — V [-e^^tanhr]'" Inn) 



6, out 



n=0 



(103) 



^e"|nn) 



6, out -I 



n=0 



where, in our case, = 7r/2 with the squeezing operator 
defined as 5 (re^^) = exp[-r(e^^at6t - e-^^a6)]. The 
squeezing parameter r can be obtained using the relation 
(6t6) = (coshr-l)/2, as 



= arcosh ( 2e^^'^/'" - 1 



(104) 



1/(1 



The t elep ortation fidelity is given by F 
e-2-) H. 

Let us discuss the fact that the entangled state in the 
original frame is given by I^6(/^)^out(<^out)|TMS)5,out. 
Here l^out(<^out) is the displacement operator of the out- 
put mode, which is displaced by o^out as a consequence 
of the displacement of the output operators a^^iu) by ol^^ 
(analogously to the discussion in App. ID 2^ . First, let us 
generally define the teleportation scheme as the map A, 
such that it teleports a light state \i))e as follows 



A [|TMS)6,out O |V')e] = 



(105) 



Here, the subindex e labels the external system contain- 
ing the state that one wants to teleport. Let us re- 
mark that perfect teleportation = 1 can only 
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be achieved for the maximahy entangled state r ^ oc. 
In order to determine the output state in the original 
frame, let us first transform the initial state 



CO 



6, out 5 



n=0 



(106) 



where O5 = E^^o The relation A 

B'^^lnn) = AB'^ (g) H-Xlnl^^) been used, where 
B'^ denotes the transpose of B. Using this relation, the 
output state of the teleportation scheme with the original 
state is given by 

A [I^6l^out|TMS)5,out ^ \^)e] = DbObDl^,\^)i,. (107) 

This gives the final state of the teleportation protocol in 

the original frame. Note that I^Q^^(aout) = ^lut(^out)- 
Therefore, one can get rid of this displacement by tele- 
porting the state D{a'^^^)\ilj)e^ such that the state tele- 
ported in the mechanical system is given by I^^(/3)0|?/;)5 
(the displacement Di){/3) can also be reduced by vary- 
ing the center of the trap when switching off the cav- 
ity lasers). Besides, note that one can in principle also 
choose the appropriate initial state in order to pre- 
pare a desired mechanical system |(/)), such that Di)\(j))i) = 
DM),. 



VI. MECHANICAL TOMOGRAPHY BY TIME 
OF FLIGHT 



After showing how to prepare non- Gaussian states, in 
this section we discuss how to measure them. In gen- 
eral quantum optomechanical systems, the proposal is to 
transfer the mechanical state into a well controlled quan- 
tum system, such as a qubit, and probe the state in that 
system. This could be analogously done in our setup by 
mapping the mechanical state to the cavity mode by us- 
ing the enhanced beam-splitter interaction and perform 
full tomography of the output field. However, this tech- 
nique would suffer from the drawback that the output 
field would contain a quantum state displaced by the 
large driving field and therefore, the signal-to- noise ra- 
tio would be challenging for experimental detection with 
present day technology. 

In this section we propose an alternative method to di- 
rectly perform full tomography of the mechanical system. 
In particular, we exploit the analogy of levitated nanodi- 
electric objects to atomic physics, more specifically to 
cold gases, where time of fiight measurements are used 
to experimentally probe different many-body states [6|. 
The time of fiight protocol to perform direct full tomog- 
raphy of the mechanical state is the following (see Fig.[6|): 





z{t + tf) « {tf - t)p{U)/M 



FIG. 6: Schematic illustration of the time of flight proto- 
col to perform full tomography of the mechanical state. The 
momentum operator at times te, which corresponds to the 
rotated phase quadrature xi^tte + 7r/2), is determined by 
measuring the position of the dielectric after some time of 
flight. By repeating the experiment at different times te, one 
can perform full tomography of the mechanical state. 



1. We consider that at t = a particular state 

in the mechanical system is prepared (for instance, 
a non- Gaussian state using the light-mechanics in- 
terface introduced in Sec.[Vj). Just after the prepa- 
ration of the state, the cavity field is switched off 
and only the optical trapping remains on. During 
these transient times the center of the trap has to 
be changed in order to account for the variation in 
the classical force created by the driving field, as 
discussed in the light-mechanics interface. 

2. Then, during some given time tg, the system is 
evolving within the harmonic potential, such that 
the mechanical momentum operator in the Heisen- 
berg picture is given by 

= ip^(6te-«* - fte-'"^'*), te[0,te], (108) 

where Pm = (Mwt/2)i/2. 

3. At t = te, the trap is switched off and the nan- 
odielectric falls freely during the time of flight tf, 
such that the distance from the center of the cavity 
along the cavity axis is given by 



Z{te + tf) = Z{te) + {tf - te) 



Pjte) 
M 



pjte) 
M ' 



(109) 



where we assume that is sufficiently large such 

that {tf - te)p{te)/M > Z{te). 

4. At t = te -\- tf^ the position z{te -\- tf) is measured 
{e.g. by imaging the object and measuring the cen- 
ter of the light spot in the screen), which means 
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that the in-trap momentum p{te) is effectively mea- 
sured. 

5. The experiment is repeated in order to obtain 
statistics for any time te G [0,27r/a;t]. 

The key observation is that with the data obtained in 
this protocol, one has the statistical distribution of the 
rotated quadrature phase operator 



(110) 



which permits the reconstruction of the Wigner func- 
tion [iHI, [IQ| , and therefore contains all the information 
about the mechanical state |^) (of course, it does not 
have to be a pure state). Indeed, there exists the follow- 
ing one-to-one relation between the momentum operator 
and the rotated quadrature phase operator. 



(Ill) 



Let us now discuss some experimental considerations. 
First, we will estimate the order of magnitude of (and 
therefore the time of flight distance df = where 
g is the gravitational acceleration). In particular, let us 
assume that after the time of flight the position can be 
measured with a resolution given by Sz. This implies that 
the object has to spread over a distance much larger than 
6z^ which means that tf M6z/pm is required. Using 
the parameters given in App. [El one obtains that tf is 
of the order of tens of ms, which would require a time 
of flight distance of the order of one cm. Although this 
position resolution is very feasible, the requirement could 
be even relaxed with the same duration of time of flight. 
The idea is to amplify the oscillation via driving the field 
with a blue-detuned laser prior to letting the object fall. 
More specifically, let us assume that just after the prepa- 
ration of the mechanical state, one impinges the cavity 
with a laser detuned to the blue sideband of the cavity. 
This corresponds to including an additional step (point 
l.b) between steps 1 and 2 in the previous protocol. The 
blue-detuned driving is performed during a certain time 
T < 1/r (where F is the decoherence rate when the cav- 
ity field is switched on). After this amplification, the 
momentum operator is transformed into 

p{r) = ip^p+(r)(6te-*- - be-'^n + Pcav(r), (112) 

where p+ (r) is the amplifying parameter given by 



-Kt/2 



cosh(xr) 



2X 



sinh(xr) 



(113) 



with X = v^+^V^- The term Pcav(T) results from 
the entanglement of the mechanical system to the cav- 
ity field due to the two mode squeezing interaction. It 
reads Pcav{r) = {q'' {T)e'^^ai{0) + H.c), where q{t) = 
— i^e~'^^/^ sinh(x^)/X7 fulfills (Pcav('?")) = (the cav- 
ity field is empty at t = 0) and {pl^^{r)) = \q{r)\'^. 



After this amplification, step 2 of the protocol follows. 
If one assumes g = n = 2ti x 100 kHz, and r = 0.02 
ms, one obtains that P+(t) ~ 10^ and hence with the 
same time of flight tf the required resolution is only 
5z <C tfPmP+{T) /M 100/im; three orders of magni- 
tude lower. Note that the amplification is restricted by 
keeping the nanodielectric object in the region, where it 
still sees the slope of the standing wave, i.e., the condition 
XoP+(r) < 1 nm has to be fulfilled, where Xq ~ 10~^^ m 
is the ground state size. 

Let us remark that the rotated quadrature x(6^) could, in 
principle, also be measured by a quantum non-demolition 
measurement. This could be done by using the back- 
action evasion scheme proposed by Braginsky in the 80 's 
[7l|, and recently revised from a quantum noise per- 
spective [72|. This protocol would also benefit from the 
prominent property of levitating objects; being free of 
any thermal contact. The key idea of this method is 
to impinge the cavity at the two motional sidebands, 
a scheme that has already been realized with trapped 
ions [73|,[7|. 

The time of flight protocol presented in this section ex- 
ploits the unique property of using levitating objects in 
quantum optomechanical systems; the mechanical res- 
onator is unattached to other objects and therefore can 
fan. 



VII. CONCLUSIONS OVERLOOK 



We conclude by summarizing and giving an overlook of 
the contents presented in this article. First, we have de- 
veloped a quantum theory to describe the coupling of 
light to the motion of dielectric objects inside a high fi- 
nesse optical cavity. The main result is the derivation 
of a master equation describing the joint state of the 
center of mass motion and the cavity field. In paral- 
lel, we have derived a quantum elasticity theory to show 
that the center of mass decouples from the internal vibra- 
tional modes for sufficiently small objects. This theory 
has been applied to describe the experimental proposal of 
using an optically levitating nanodielectric as a cavity op- 
tomechanical system [l|, Q • The master equation allows 
us to describe the coherent dynamics as well as the dissi- 
pative processes. More specifically, we have obtained the 
decoherence rate for the mechanical mode, the enhanced 
cavity decay rate due to light scattering, and a renormal- 
ization of the light-matter interaction Hamiltonian due to 
virtual photon exchange processes. This theory supports 
the statement that the center of mass motion of a levi- 
tating nanodielectric inside an optical cavity behaves as 
a mechanical resonator of very high quality. 

In the second part of the article, we have developed a 
light-mechanics interface with the aim of bringing lev- 
itating objects into the quantum regime. This can be 
used to prepare non-Gaussian states such as superposi- 



21 



tions of Fock states. We have provided three protocols 
with different properties; first, the reflected one-photon, 
which requires strong coupling and a measurement of 
the output field. Second, the perfect mapping, which 
does not require the measurement by time-modulating 
the optomechanical coupling. Third, the teleportation 
in the bad cavity limit, which can be used in the weak 
coupling regime. These light-mechanics interfaces apply 
to other optomechanical systems and provide an effec- 
tive way to obtain non-linearities. Besides, these input- 
output protocols required a formalism to be described in 
the Schrodinger picture in order to obtain the final state. 
Finally, we have proposed a method to perform direct 
full tomography of the mechanical state. This method 
exploits the levitation of the mechanical resonator, since 
it consists in measuring the center of mass position after 
letting the object fall. The position after time of flight 
carries information about the momentum in the trap, 
and by repeating the measurement at different times of 
evolution in the harmonic trap, one can perform full to- 
mography. 

The theory and the protocols introduced in this article 
apply to a large variety of setups and dielectric objects, 
even to microorganisms [l^. This work opens many fur- 
ther directions that we are currently investigating. Some 
of them are the following: (i) the possibility to apply time 
of flight experiments to prepare and measure macroscopic 
superpositions of the levitating object, that is, states in 
which the object is in two macroscopically distant (larger 
than its radius) positions, (ii) A thorough study of higher 
order light scattering processes in dielectric objects, (iii) 
The possibility to circumvent the decoherence processes 
due to light scattering by using magnetic levitation of 
micron objects, (iv) To reduce light scattering by using 
dielectrics of other shapes, (v) To address the internal 
modes of the object for high frequency resonator optome- 
chanical purposes, (vi) To use objects with internal de- 
grees of freedom, such as nanocristals with NV centers 
in order to couple the internal degree of freedom to the 
center of mass mode. 

As mentioned in the introduction, the project of cavity 
optomechanics with levitating objects aims at applying 
the quantum techniques developed to control and manip- 
ulate atoms back to the nanodielectrics that were first 
used by Ashkin. It is our hope that this article stim- 
ulates further theoretical and experimental research in 
this direction. The ultimate goal of these investigations 
is to explore the boundaries of quantum mechanics, which 
may reveal unexpected, and fascinating new insights. 
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Appendix A: Light scattering equation 



In this Appendix, we show how from the light-matter in- 
teraction Hamiltonian Eq. one can derive the scat- 
tering equations which can be used to compute the total 
electric field inside the dielectric object. We decompose 
the total electric field into 



i?(x) = 



(27r)3/2 



/ 



ikx 



a(k)-H.c.), (Al) 



where uj\^ is the frequency of the different light modes and 
eo the vacuum permittivity. Starting from the Hamilto- 
nian (f24|) (including the free term Hj^^^, Eq. (fT8|) . of the 
modes a(k)), one can connect the electric field E{x.) to 
the electric field without the presence of the dielectric ob- 
ject, Eo{-k). To achieve this, we determine the equation 
of motion for the annihilation operator a(k) 

d{k,t) = -icka(k,t)-^^^yg J d'xE{^,t)e''^-. 

(A2) 

In order to obtain £^(x), one needs to formally integrate 
Eq. (|A2|) over time, multiply both sides of the equation by 
i\/^/c/2eo/(27r)^/^e~^*^^, and integrate it over k. Given 
that the electric field is varying with the laser frequency, 
one can move to the rotating frame, where £^(x, t) is 
slowly varying, E^((^^t) = e^^'^^^£^^(x). This justifies 
the assumption E^{-x^t) ^ E^^^^t') which permits to 
simplify the integration over time, leading to 

^+(x,t) =^o^(x,t) 

J J 2(27 



^+(x' 



Jo 



2(27r)3 

^^g-i(u;z.-u;,)T 



(A3) 



•^-(x^^) / dre^^'^^+'^^)^ 
Jo 



The terms containing E~ (x, t) are rotating fast compared 
to any other time scale and will hence be neglected within 
a rotating wave approximation (RWA). Carrying out the 
integration d^k in Eq. (\A3\i gives us a function that de- 
cays very quickly in r. This permits us to extend the 
upper integration boundary t to oo and hence to obtain 



E{^,t) =Eo{^,t)^ap J d^x'a(x ,x)£;(x^t), 



(A4) 



where ^(x,x') = /c^ sin(/cL|x— x'|)/|x— x'|. This equation 
has the same structure as a scattering equation where 
^(x, x') takes the role of the propagator. 
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Let us remark here that in order to determine £^(x, t) 
for any size and shape of the object, the scattering wave 
equation (jA4p has to be solved. However, this equation 
can only be solved approximately (see [75|, |76| for a so- 
lution). Merely in the special case of a spherical object, 
the electric field can be determined exactly by expansion 
in spherical waves, the Mie solution [t^, [l8|. Indeed, this 
solution coincides with the discrete-dipole aproximation 
for perfectly spherical objects [76[. In the limit of very 
large spheres, ^ A, a ray optics approach has to be 
used to determine the forces on the sphere [TOf . 



Appendix B: Optomechanical parameters 



In this Appendix we show how from the light-matter in- 
teraction term of the total Hamiltonian Eq. ([24|) . one can 
easily obtain the optomechanical Hamiltonian discussed 
in Sec. ImC!! 



1. External trapping with optical tweezers 



The trapping of the sphere can be achieved either by 
using optical tweezers [49| or two optical cavity modes 
[H, [2[. For the optical tweezers, we assume a Gaussian 
beam 

[^tw(x)] = i.o^exp(-^), (Bl) 

where = [Pt/i^oCTrW^)]^^'^ , Pt is the laser power, Wt 
is the laser beam waist, W{y) = Wt[l^ {yX/ {7rW^)f]^^'^ 
and assume the beam is aligned as sketched in Fig. [H 
The interaction between the sphere and the light field is 
described by Eq. ([24|) . If the object is smaller than the 
laser waist and is placed close to the beam center, one 
obtains after integrating over a sphere of radius i?, mass 
M, density p and relative dielectric constant e^, that the 
Hamiltonian is of the form of a harmonic oscillator with 
frequency Ut given by 



4ec J_ 
pc 



pc 



(B2) 



corresponding to the x— and z— direction in our setup. 
Here, / is the field intensity, the laser waist can be ap- 
proximated by Wt ^ A/ (ttA/') , M is the numerical aper- 
ture, and k = 27r/A the wave vector. In the direction of 
light propagation, the direction in our configuration, 
the trapping frequency is 



26c 

pc ' 



(B3) 



which is reduced by a factor of A/'^/2 compared to the 
trapping in the x— and z— direction. This lower trap- 
ping frequency can be enhanced by the use of a second 



optical tweezers perpendicular to the first one. Besides, 
the scattering force will change the equilibrium position 
of the object in the direction of light propagation. To 
circumvent this, a second tweezers of the same intensity 
and waist, but with a different polarization, can be used. 



2. The optomechanical coupling 



The optomechanical coupling arises from plugging the 
cavity mode into Eq. ([24|) . For small spheres, we choose 
a TEM 00 mode as the cavity mode in the presence of 
the sphere. Given that the sphere has a radius smaller 
than the laser waist and is placed close to the center of 
the cavity, one can approximate the square of the elec- 
tromagnetic field close to the center of the beam by 



[^cav(x)]' 



1 



2(^2 



^ cos^ {kcZ — (f) a^a. 



(B4) 

Here, Vc = ttW^L/A is the cavity volume, Wc = 
[AL/(27r)]^/^ is the laser's waist at the center of a confo- 
cal cavity, L the cavity length, ujc the cavity's resonance 
frequency, A the laser wavelength and a{a^) the annihila- 
tion (creation) operator of cavity photons. Furthermore, 
we presume that the laser is aligned such that the wave 
vector of the cavity mode kc points in z-direction. The 
integration over the volume V around the center of mass 
position r = (x, z) leads to 



1/ 

^ JV{v 



with 



/(r) = 



Vc, [W^ - 2{x^ + y^)] cos^ 



(B5) 



(B6) 



where V is the volume of the sphere, z the center of mass 
position in z-direction and a phase shift. If the sphere 
is trapped at the maximum slope of the standing wave, 
that is, xo = ^0 = and = 0, = 7r/4, it is justified to 
expand Eq. (jB6p to first order in the 2:-coordinate. The 
zeroth order contribution leads to a constant shift of the 
trapping frequency given by 



(B7) 



which we will always use in the article without explic- 
itly denoting the tilde. After quantization of the z- 
coordinate, the first order contribution yields the op- 
tomechnical coupling g^a^a(\)^ + 6), where the optome- 
chanical coupling reads 



90 



^/2Mujt cVc 



(B8) 
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3. Shift due to interference of the tweezers and the 
cavity field 



In this section, we show that the term of the shift Hamil- 
tonian H^h discussed in Sec. IIIIDI yields a shift in the 
trapping frequency as well as in the equilibrium position 
of the dielectric object. We first recall the term leading 
to this shift, 

Hsh = -eceo / c^^xfcav(x)ftw(x), (B9) 

JV(r) 



where 



4av(x) = 2 J cos{kcZ - (p)\a\. 



(BIO) 



Then, by expanding up to second order around the equi- 
hbrium position of the center of mass coordinate, one 
obtains 

H,h=io,i,b^b + U{b^ +b), (Bll) 
where the parameters are given by 



2 _ ec\a\ cOcPt k,Wf + 2 



p V VcCTT wi 



Csh = -ecV\a\kcZo^ 



I COcPt 

VeCnW^ 



(B12) 



Therefore, the total trapping frequency of the object is 
given by cj^ = uJt-\-0Jsh- The change in the equilibrium po- 
sition is considered when displacing the mechanical mode 
in Sec. IVAl 



Appendix C: Ground state cooling 



Using the theory of ground state sideband cooling in op- 
tomechanical systems, p5l-[29|. the minimal number of 
phonons attainable is given by 



4wt, 



(CI) 



Taking into account heating mechanisms, such as the re- 
coil heating due to light scattering Fsc, and others Pothers 
{e.g. Brownian motion heating due to the surrounding 
gas, laser noise, blackbody radiation, etc. [H, Q), which 
can be shown to fulfill Pothers <^ Psc, the final phonon 
occupation number is given by 



" Pothe r 



(C2) 



For g < tv, the maximal achievable cooling rate is given 

by 



r_ =4 



(C3) 



where the detuning is A ~ cj^. By using the experimental 
parameters described in App. [El nu can be made much 
smaller than one. 



Appendix D: Displacement of the output modes 



In this section we show how the expression of the dis- 
placement of the output modes, a^j^ appears naturally 
by computing the steady state obtained when the driv- 
ing field is switched on. Then we discuss how to measure 
a photon created on top of the coherent cavity field in the 
output field. To simplify the problem we assume a cavity 
of resonance frequency ujc^ driven by a laser at cjl- In 
the rotating frame at the laser frequency, and by defining 
A = ujc — ujl, the Hamiltonian reads 



H =Aah 



/oo 
ujal 



{uj)ao{uj)duj 



7(cj)(a^ao(cj) — Il.c.)duj. 



1. Steady-state with a driving field 



The initial state of the system is assumed to be 



(Dl) 



m 



a)® / S{uj)D{ao)\n), 



(D2) 



that is, the cavity state is in a coherent state with phase 
a, and all the output modes are empty, only the laser 
mode is in a coherent state with phase ao (which is re- 
lated to the laser power). In the following, we aim at 
computing the final state |st) = lim^^oo exp[— iJYt] |in). 

First, let us write the Heisenberg evolution equations for 
the cavity mode and the output modes: 



/oo 
7(cj)ao(cj, t), 

do{uj,t) = —iujao{uj^t) — ^{uj)a{t). 



(D3) 



Then, one can formally integrate the differential equation 
for ao(cc;, t) 



^ao(cj,0) - 7(0; 



) / dTa{T)e 
Jo 



0{t-T) 



(D4) 

This solution can be introduced into the differential equa- 
tion for a{t). By using the approximation jioj) ^ 7(0) = 
Y^a^/tt, one gets 



d{t) = -(iA- 



J — c 



j{uj)e-'^'ao{co,0), (D5) 
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which can be triviahy integrated to 
a{t) =e-(^^+")^a(0) 

/oo 
7(w)e-''^^ao(w,0)e-(''^+'')(*-^). 

(D6) 

By taking the average value of this expression, using that 
(a(0)) = a and (ao(cj,0)) = ao, one gets 

1 _ p-iiA+K)t 

{a{t)) = e-('^+'^)*a + ^{0)ao^^ . (D7) 

In the steady state, one obtains 

7(0)ao 



a = hm (a(t)) 

t^oo^ ^ ^' iA + iA + 



(D8) 



Note that we have assumed that the initial coherent state 
of the cavity is equal to the steady state obtained when 
driving the cavity with the laser. We have also related 
Qfo to the usual frequency VLl = y^^J^J^^ Pc being the 
laser power. Let us now compute the mean value of the 
output modes, which after some algebra is given by 

{ao{w,t))=ao5{u:)-^{w) f dT{a{T))e-'^^'--^ 

(D9) 

Jo 

Then, the steady state phase of the output modes can be 
expressed by 



Qfo; = lim (ao(cj,t)) 

t^oo 



(DIO) 



which coincides with the expression used in Eq. (|6Q]) . 

It is trivial to show that the Hamiltonian is invariant 
under the displacement operation V = Z^^I^out, such that 
DlaDa = a + a, and Dl^^ao{uj)D out = ao{uj) + o^c^. By 
usmg that V uj-^duo = 0, one can check that 

V'^HV = H. (Dll) 
This implies that the steady state 

/CO 
dujD{au;)\On) (D12) 

is indeed an eigenstate of the Hamiltonian: 

i^|in) = VV^HV\m) = VH\On) = 0. (D13) 

2. Measurement of a photon 



In this section we want to compute the displacement of 
the output mode of the cavity. To do so, we assume 



that at t = a photon is present inside the cavity in the 
displaced frame, that is. 



(D14) 



Therefore, using the Wigner-Weisskopf formalism, one 
can obtain the state at some later time 

nOO 

\ilj{t)) = Ca{t)\in) ^ / dujc{uj,t)al{uj)\On), (D15) 



where the coefficients are given by 



c{uj^ t) 



(D16) 



i{uj - Ao) + 



For large t, the final state is given by \i^{t)) = aI^^ tl^^)' 
where the collective output mode is defined as Aout,t = 
/ (l)out{^)e^^* ao{uj)duj ^ with the mode function 



0out(^) 



7(cj) 



Hi — — Ao) 



(D17) 



Let us now compute how many photons will be encoun- 
tered in this collective photon mode after transforming 
back to the non-displaced frame. By using the expres- 
sion of the displacement of the output modes a^^ one 
can obtain after some careful manipulation that the dis- 
placement of the output mode Aout,t7 defined as 



(D18) 



/OO 

is just given by c^out = 

3. Switching off the driving field 



In this Appendix we want to discuss the final state of 
the one-photon protocol once the driving field has been 
switched off. The Hamiltonian in the frame rotating with 
the laser frequency cj^,, is given by 

/oo 
uja\{ijj)aQ{uj)duj 
-OJT. 



a(6'^ + 6) 



A(t)(6t+6), 



'y{uj){a^ ao{uj) — H.c.)^^^^ 



(D19) 



where the term with X{t) accounts for the variation of 
the center of the harmonic trap. By writing the Langevin 
equations, and considering that there is no input fields 
since they have already been switched off, one obtains 



d{t) = — iAa(t) — K:a{t) — igoa{b^ + 6), 
b{t) = —iuotb — igoa^a — iX{t). 



(D20) 
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By displacing the operators by a' = a+a(t), and choosing 
the restriction 



a = — iAa — K,a^ 
= -igo\af -iX{t), 

one obtains the fohowing equations: 

d{t) = -iAa{t) - Ka{t) - i^o|<^|(^''' + 
h(t) = -\uJth - \gQ\a\{o) + a). 



(D21) ^ 



(D22) 



In the interaction picture one can perform the RWA in 
order to get 



a(t) = —na(t) — ig{t)b, 
h{t) = -ig{t)a, 

where 

g{t)=g^\ame-^'. 
The equation for h(t) is then given by 



(D23) 



(D24) 



(D25) 



This can be solved in order to predict the variation of 
mechanical state by switching off the driving field. 



Appendix E: Experimental parameters 



In order to illustrate the experimental feasibility of the 
proposal, we will choose a set of experimental parame- 
ters. 

1. Dielectric object: We assume nanospheres fab- 
ricated of fused silica with a radius R = lOOnm, 
density p = 2201kg/m^ and a dielectric constant 
^[cr] = 2.1 and ^[e^] - 2.5 x 10"^^ Their Young 
modulus is y = 73 GPa and their Poisson constant 
a = 0.17, giving internal vibrational modes with 
frequencies of the order of ~ 10^^ Hz. 

2. Cavity: We assume a confocal high finesse cavity 
of length L = 4 mm and finesse = 5x10^ leading 
to a cavity decay rate = C7r/2TL = 27r x 44kHz. 

3. Lasers: The optical tweezers is constructed with 
a laser of power Pt = 15mW at a wavelength A = 
1064 nm and a lense of high numerical aperture 
JV = 0.8. The cavity is impinged by a laser of 
power Pc = 0.1 mW, wavele ngth A = 1064 nm, 
which gives a waist of Wc = ^/Xd/27T ^ 26 /im. 

4. Dissipation due light scattering: The decay 
rate of the cavity due to the light scattering is tZsc = 
27r X 6 kHz (see Fig. [3 for the dependence with 
the radius of the sphere). The mechanical motion 
decoherence rate is given by Fgc = 27r x 9 kHz. 
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FIG. 7: The ration between the decay rate of the cavity due to 
light scattering Ksc over the standard decay rate of a cavity of 
finesse = 5 x 10^ is plotted as a function of the radius of the 
sphere. As can be observed, photon losses due to scattering 
are dominating for objects larger than 150 nm. 
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FIG. 8: Final phonon occupation number using spheres of 
radius R and for the parameters given in the text. Ground 
state cooling is possible for spheres smaller than R ^ 200 nm. 



5. Optomechanical parameters: The tweezers 
supplies a harmonic trap for the object of frequency 
ujt = 27r X 136 kHz in the transversal direction 
and cj|| = 27r X 77 kHz in the direction of light 
propagation. The steady state photon number is 
|ap ~ 3.7 X 10^. The optomechanical coupling to 
the CM degree of freedom of the sphere is given by 
^0 ^ — 27r X 7 Hz, which is enhanced by a factor 
of \a\ to g = — 27r x 131 kHz. The frequency of 
the cavity photons is given by c^c = 27r x 2.8 x 10^^ 
Hz. Note that uJt/{i^-\- i^sc) ^ 3 (condition for the 
RWA), g / {tz -\- tisc) ^ 1 (the strong coupling regime 
used in the light-mechanics interface), g/T^c ^15, 
and uJt/T sc ~ 17 (number of coherent oscillations). 

6. Ground state cooling: The cooling rate is given 
by F_ = 27r X 1.6 MHz. The final occupation num- 
ber is given by um = 0.014, see Fig. [8] for the de- 
pendence with the radius of the sphere. 
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